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are reviewed. The impact that a direct observation of this exotic process would 
have on elementary particle physics, nuclear physics, astrophysics and cosmology 
is profound. Now that neutrinos are known to have mass and experiments are be- 
coming more sensitive, even the non-observation of neutrinoless double-beta decay 
will be useful. If the process is actually observed, we will immediately learn much 
about the neutrino. The status and discovery potential of proposed experiments 
are reviewed in this context, with significant emphasis on proposals favored by re- 
cent panel reviews. The importance of and challenges in the calculation of nuclear 
matrix elements that govern the decay are considered in detail. The increasing sen- 
sitivity of experiments and improvements in nuclear theory make the future exciting 
for this field at the interface of nuclear and particle physics. 
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I. INTRODUCTION 

Neutrinoless double beta decay (/3/?(0i/)) is a very slow lepton-number-violating nuclear 
transition that occurs if neutrinos have mass (which they do) and are their own antiparticles. 
An initial nucleus {Z,A), with proton number Z and total nucleon number A decays to 
{Z + 2, A), emitting two electrons in the process. A related transition, called two-neutrino 
double-beta decay (/3/5(2z/)) results in the emission of two electron antineutrinos in addition 
to the electrons, and occurs whether or not neutrinos are their own antiparticles. (3j3{2v) has 
in fact been observed in a number of experiments. With the exception of one unconfirmed 
observation, on the other hand, (3(3{Qv) has never been seen, and searches for it are ongoing 
in a number of laboratories around the world. Other even slower and more exotic process, 
including double-positron decay, double-electron capture, and neutrinoless decay with the 
emission of a hypothetical particle called the Majoron (/?/?(0z/, x)), have likewise never been 
observed. 
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The development of effective-field theory and grand-unification schemes in the late 1970's 
and early 1980's led to the expectation that neutrinos are identical with their antiparticles 
and have nonzero mass, and engendered renewed interest in f3f3{0i>) experiments. More 
recently, neutrino-oscillation experiments have yielded compelling evidence that the three 
observed flavors of neutrinos are not mass eigenstates but rather linear combinations of those 
eigenstates (at least two of which have nonzero mass eigenvalues). These experiments also 
allow the electron neutrino to mix significantly with the heaviest of the mass eigenstates. If 
it does, the effective neutrino mass will be large enough that /3/3(0z/) may well be observed in 
experiments currently under construction or development. An observation would establish 
that neutrinos are "Majorana" particles = V, roughly speaking), and a measurement of 
the decay rate, when combined with neutrino-oscillation data, would yield insight into all 
three neutrino-mass eigenstates. This article is motivated by the recent developments in 
neutrino physics and by the interest among physicists throughout the world in a coherent 
experimental /3/3(0z/) program. 



A. The Early History 



Double-beta decay was first considered in a 1935 paper by Maria Goeppert-Mayer 



( iGoeppert-Mayeii . Il935l ). The author, who acknowledged Eugene Wigner for suggesting the 
problem, derived an expression for the decay rate and estimated a half-life of ~ 10^^ y for a 
decay with the emission of two electrons and two anti- neutrinos (l^), carrying about 10 MeV 



of energy. Two years later Ettore Ma, 



orana formulate d a theory of neutrinos in which there 



was no distinction between u and u (jMajoranal . Il937l ). and suggested an experimental test 
of his hypothesis with a reaction similar to Vj^+f^C I — Ar + e~, which was later searched 



for (and not found) by Raymond Davis (iDavis Jr.l . Il955l ). I t was Giulio Racah, however, 
who first suggested testing Majorana's t heory with Pp (Oi') (IRacahl . 119371 ). In 1939 Furry 



calcu lated approximate rates for /3/3(0z/) (IFurryl . Il939l ). and in 1952 Primakoff (jPrimakofl, 
19521 ) calculated the electron-electron angular correlations and electron energy spectra for 
both I3[3(2v) and /5/3(Oz/), producing a useful tool for distinguishing between the two pro- 
cesses. These early contributions set the stage for many years of experimental and theoretical 
activity. 



The review by iHaxton and Stephenson Jr.l (1l984j ) contains a chronology of experiments 
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from 1948 through 1983. There were some early claims of observation. iFiremanI (119491 ) 
repor ted observing t he (3(3 of ^^^Sn in a laboratory experiment, but retracted the claim 



later (IFireman 



195^ . The first geochemical observatio n of with an estirnated h alf-life 



1.4 X 10 y, was reported in 1950 (llngram and Reynolds 



195q[) 



The 



of Tf/2 (''°Te) 

first actual laboratory observation of (3(3{2v) was not made until 1987 (lEUiott et ali 119871 ). 
Here, we concentrate on experiments developed since the late 1980s, referencing earlier work 



where a ppropriate. The early dey e 



example 



Avignone and Brodzinski ( 



mm-. Elliott and Engell f2004h: 



Haxton and Stephenson Jr. 



Suhonen and Civitaresd (119981 ): 



Zdesenk 



3 (2q3) 



opme n ts have been co v ered well in other reviews, for 



1988): 



Avignone al. 



1984): 



Elhott and Vogell (12002 



Moe and Vogell (11994h: 



(2005); 



Doi et al. 



1985 



aessler and Simkovic ( 



Primakoff and Rosen 



Tomodal (Il99lh : iTretvak and Yu G. Zdesenkol mm 



Eiiri 



1998) 



1981) 



B. Overview of Theory and Recent Experimental Developments 

A typical (3(3 candidate is an even-even nucleus (Z, A) which pairing forces make more 
bound than its (Z+1,A) neighbor, but less so than the (Z+2,A) nuclide, as shown in Fig. 
[U In Fig. [2] we depict (3(3{2v) and neutrino-exchange-driven (3(3{Qv). The rate of (3(3{2v), 
which has been measured in 10 isotopes (see Table [IT|) . can be written as 

{Tlj^)-' = GUQpp,Z)\M^^\ (1) 

where G2u{Qi3f3, Z) is the four-particle phase space factor, and is a "nuclear matrix 
element" for this second-order process. This decay conserves lepton number, does not dis- 
criminate between Majorana and Dirac neutrinos, and does not depend significantly on the 
masses of the neutrinos. The rate of /3/3(Oz/), if driven by the exchange of light Majorana 
neutrinos, is approximately 

(Tf/^)-' = GMpp,Z)\M^,Wmppf , (2) 

where GQ^iQpp, Z) is the phase space factor for the emission of the two electrons, Mq^ is 
another nuclear matrix element, and {mpp) is the "effective" Majorana mass of the electron 
neutrino: 



5 



k 



(3) 



Here the m^'s are the masses of the three hght neutrinos and U is the matrix that transforms 
states with well-defined mass into states with well-defined flavor [e.g., electron, mu, tau). 
Equation [2] gives the /3/3(0z/) rate if the exchange of light Majorana neutrinos with left handed 
interactions is responsible. Other mechanisms are possible (see Sections IIIII and lIV.DI) . but 
they require the existence of new particles and/or interactions in addition to requiring that 
neutrinos be Majorana particles. Light-neutrino exchange is therefore, in some sense, the 
"minima" mechanism and the most commonly considered. 

That neutrinos mix and have mass is now accepted wisdom. Oscillation experiments 
constrain U fairly well — Table [T] summarizes our current knowledge — but they determine 
only the differences between the squares of the masses rrik {e.g., m\ — m\) rather than the 
masses themselves. It will turn out that /?/5(0z/) is among the best ways of getting at the 
masses (along with cosmology and /3-decay measurements), and the only practical way to 
establish that neutrinos are Majorana particles. 

To extract the effective mass from a measurement, it is customary to define a nuclear 
structure factor Fjq = GQ^iQ^p, Z)|Moi.pmg, where rrie is the electron mass. (The quantity 
Fm is sometimes written as Cmm-) The effective mass {rripis) can be written in terms of the 
calculated Fn and the measured half life as 



(m^^)=me[F^T0;2]-i/2. (4) 

The range of mixing matrix values given below in Table [H combined with calculated values 
for F/v, allow us to estimate the half-life a given experiment must be able to measure in order 
to be sensitive to a particular value of {mpjs). Published values of Fn are typically between 
10"^'^ and 10~^^ y~^. To reach a sensitivity of {rriisf^)^ 0.1 eV, therefore, an experiment must 
be able to observe a half life of 10^^ — 10^"^ y. As we discuss later, at this level of sensitivity 
an experiment can draw important conclusions whether or not the decay is observed. 



The most sensitive 



V fjAalseth et al. 



(lArnaboldi et al. 



imits thus far ar e from the Heidelberg-Moscow experiment: 



Baudis et al. 



19991), the IGEX experiment: T[/^ 



•Ov /76 



Ge) > 1.6x10 



25 



2004J) . and the CUORICINO experiment T^J^i^^^Te) > 3.0 x 10^^ y 



2002a! J 

2OO5I, 2OO7I). These experiments contained 5 to 10 kg of the parent isotope 
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and ran for several years. Hence, increasing the half-life sensitivity by a factor of about 100, 
the goal of the next generation of experiments, will require hundreds of kg of parent isotope 
and a significant decrease in background beyond the present state of the art (roughly 0.1 
counts/ (keV kg y). 



It is straightforward to derive an approximate analytical expre ssion for the ha 



which an experiment with a given level of background is sensitive (lAvignone et al. 



f-hfe to 



20051): 



416x_^/eaX / Mt 



Here n^- is the number of standard deviations corresponding to a given confidence level 
(C.L.) — a CL of 99.73% corresponds to no- = 3 — the quantity e is the event-detection and 
identification efficiency, a is the isotopic abundance, W is the molecular weight of the source 
material, and M is the total mass of the source. The instrumental spectral- width A{E), 
defining the signal region, is related to the energy resolution at the energy of the expected 
PPlOu) peak, and b is the specific background rate in counts/(keV kg y), where the mass is 
that of the source, as opposed to the isotope. Equation ([5]) is valid only if the background 
level is large enough so that the uncertainty is proportional to \JhK{E). For a 200-kg ^^Ge 
experiment with a background rate of 0.01 counts/(keV kg y) and an energy resolution 
of 3.5 keV, running for 5 years, the values for these parameters are Mt = lO^kg ■ y,e = 
0.95, a = 0.86, W = 76, and A{E)= 3.5 keV. This results in a Aa half-life sensitivity of 
^1/2 (4*^) "^^Ge) = 1.9 X 10^^ y. The background rate quoted above is conservative for a Ge 
experiment, only a factor of 6 below that of the Heidelberg- Moscow and IGEX experiments. 
A background rate of 0.001 counts/(keV kg y) would allow a 4cr discovery with = 
6 X 10^^ y. But an experiment with even modestly lower efficiency or poorer resolution must 
attain much lower background rates to have comparable sensitivity. 

These numbers characterize the level future experiments will have to reach to make a 
significant contribution to the field. Later, we will discuss a number of proposed experiments 
and attempt to estimate their sensitivity. 



C. The Claimed Observation 



In 2 001, a subset of the Heidelberg- Moscow collabo r ation 



In 2001, a 


2001a 




2003 





Klapdor-Kleingrothaus and Krivosheina 



2006 



( Klapdor-Kleingrothaus et al 



Klapdor-Kleingrothaus et al 



2001bl . I2OO3I ) claimed to observe evidence for a (3/3{0i') peak in the spectrum of their Ge 



experiment at 2039 keV. This claim and later pape r s by th e same group e. 



of crit i cal replies, for e x amp 



(I2QO1I); 



Zdesenko et al. 



e Refs. 



Aalseth et al. 



fl2002bh : iFerudio et al. 



icited a number 



(120021 ): 



Harney 



( 120021 ) ■ But whether or not the result is valid, the experiment was 
the most sensitive so far. The parameter values were Mt = 71.7 kg y, 6 = 0.11 counts/(keV 
kg j), e = 0.95, a = 0.86, W = 76, a nd A(E) = 3.27 keV. The nu r nber o f counts under 



the peak at 2039 keV was 28.75 ± 6.86 (iKlapdor-Kleingrothaus et ail l2004al ). Substitution 
into Eq. ([5]) yields T°/^2 Ge) = 1.6 x 10^^ y, a lifetime comparable to the claimed most 
probable value, 2.23 x 10^^ y. At least nominally, the experiment had a 4a discovery potential, 
and cannot be dismissed out of hand. Since this analysis does not account for statistical 
fluctuations, background structure, or systematic uncertainties, the actual confldence level 
could be signiflcantly different. But the only certain way to conflrm or refute the claim is 
with additional experimentation, preferably in ''^Ge. 



To this end, the GERDA experi ment is under c onstru ction in the LNGS ( lAbt et al. 



20041 ) and the Majorana project (IGaitskell et all. 120031) is 



The CUORICINO experiment in the LNGS (!Arnaboldi et all 



jeing developed in the U.S. 
200M) uses ^^°Te, and is the 



most sensitive experiment currently operating, with a lower bound of Ty2(}^^Te) > 3 x lO^"' 
y. This limit is at or near the sensitivity needed to test the 2001 claim, but uncertainty in 
the calculated value of the nuclear matrix element Mq^ (or, equivalently, Fn) will preclude 
a deflnitive statement. 

Foiled by the nuclear matrix elements! One can see even in this brief overview how nice 
it would be to have accurate matrix elements. We will address the issue of how well they 
can be calculated later. 



II. MAJORANA NEUTRINOS 

As we will see, /9/9(0z/) cannot occur unless neutrinos are Majorana particles, i.e. their 
own antiparticles. We therefore briefly review the nature of neutral fer mions. Much of what 



we say here is covered in a similar fashion but much greater depth in 



( 11987h 



Bilenky and Petcov 



We can project an arbitrary 4-spinor onto states of deflnite chirality, which are singlets 
under one of the two SU{2) algebras that make up the Lorentz algebra 5*0(3, 1). We deflne 
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the left- and right-handed projections as '^l,r = [(1 =F75)/2]^I/. Because of the minus sign in 
the Minkowski metric and the resulting need to work with \E' = \E'^7o rather than \E'^ alone, 
a Lorentz scalar cannot be constructed by contracting two-left handed spinors in the usual 
way. Instead one must contract a left-handed spinor with a right-handed one. The (scalar) 
term in the Lagranginan obtained in this way is called the "Dirac mass" term: 

jCd = -mz)^^ = -mo {^^r + W^^l) . (6) 

The terms above contract \1'2 with \l'/j (and vice versa), with 70 flipping the chirality so that 
the contraction can be made. 

If charge conservation isn't required, one can form a scalar by combining with itself 
rather than with \E'*. Such a term cannot exist in the Lagrangian for electrons or quarks 
because it destroys or creates two particles of the same charge (or destroys a particle while 
creating an antiparticle of opposite charge), but nothing forbids it in the Lagrangian of 
neutral fermions such as neutrinos. It violates lepton number, a global U{1) symmetry that 
tracks the difference between particle number and antiparticle number, but there is nothing 
sacred about global symmetries. Indeed, PP^Ou) can't occur unless lepton-number symmetry 
is violated. 

To construct a Lorentz scalar from two \[''s, we note that the "charge conjugate" of \E', 
defined up to an arbitrary phase as \E'^ = 7^\E', transforms in the correct way; its chirality is 
opposite that of \E' itself because 7^ and 75 anticommute. Thus {"^lY is right-handed, and 
we can construct left- and right-handed "Majorana" mass terms of the form 

Cm = {WlT^l + h.c) - ^ ((^rT^r + h.c) . (7) 

Although we have used all four Dirac components in this equation, it is possible that only 
the two in \E'i actually exist if there is no Dirac mass term. 

Equations (j6]) and (j7 l) can be generahzed to include N flavors. If, following 



Bilenky and Petcovi ( 119871 ) and letting be the neutrino field u, we define 

UL = , (8) 



where 
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iy^Rr 

{^rRr 



(9) 



then we find for tlie mass part of the Lagrangian, 



1- 



-Cd+m = -^{nLYMriL + h.c. , M 



Ml Ml 



(10) 



Md Mr 

The (generically complex symmetric) N x N matrices Ml and Mr are the Majorana 
mass terms and the N x N matrix Md is, the Dirac term. Except in special cases, the 
eigenstates of Cd+m will correspond to Majorana particles. The Lagrangian is invariant 
under CPT, a transformation that takes a particle into an antiparticle with spin direction 
reversed. Since the eigenstates will in general consist of 2N nondegenerate spinors, the 
components in each spinor must transform into each under under CPT, rather than into an 
orthogonal spinor, so that the neutrinos will be their own antiparticles. To see this more 
precisely, we note that M, if nondegenerate, can be rewritten in the form M = iU'^)'^ rhLl\ 
where m is a diagonal matrix with 2N positive entries irik and lA is unitary. The mass 
Lagrangian then takes the form: 



2N 



2N 



D+M 



(11) 



k=l 



k=l 



where Hl = Un'r and 



^k = <^ + {n'.^r = <Pl , (12) 

with only the cross terms surviving when the 0^ are written out as in Eq. (1121) . Clearly, 
then, the physical eigenstates (pk are Majorana particles. The interacting left-handed flavor 
eigenstates are linear combinations of the left-handed parts of these Majorana eigenstates 
(which don't have well-defined chirality). 

Dirac neutrinos are a special case, produced if, e.g., Ml and Mr are zero. We can see 
this most easily in the case of the one-flavor mass matrix where Md is simply m; 
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M 




(13) 



The eigenvalues m and -m are obtained by diagonalizing this real matrix with a real orthog- 
onal matrix. To get the two states with positive mass m we cannot use this real matrix as 
U but instead must incorporate a phase that changes the sign of the second eigenvalue: 



(14) 



The phase in this matrix that mixes the degenerate Majorana states (pi and (j)2 will cause 
the two to cancel in the neutrino-exchange diagram via which Majorana neutrinos directly 
mediate /3/3(0z/) decay. And since they are degenerate, orthogonal linear combinations xi = 
1/ \/2{(j)i + i(j)2) and X2 = l/v^(0i — "^02), that go into one another under charge conjugation, 
are also eigensates yielding Cd = (— '^/2)(0i0i + 0202) = — "^(XiXi + X2X2)- The x's make 
up the lepton-number conserving Dirac neutrino and antineutrino. 

We can also use the one-flavor case to illustrate the "see-s aw" mechanism, a nat 



ural explanation for the existeii c e of light Majo rana neutrinos (iGell-Mann et al 



Mohapatra and Senjanovic 



1980 



1979 : 



Yanagidal . Il979l ). Majorana mass terms for left-handed 



neutrinos cannot be accommodated in the standard model because those terms have dif- 
ferent quantum numbers under SU{2)l x SU{2)fj than do the Dirac mass terms. But by 
introducing new physics (e.g. new Higgs bosons with different quantum numbers) at a very 
large mass scale mn, extended models can avoid this problem. The result is often a mass 
matrix Ai (in the one-flavor example, for which Ain, Ml, and Ain become numbers mn, 
rriL, and rriu) with m^j ^ rriij ^ rriL, where m^,, which comes from our familiar Higgs vac- 
uum expectation value, is on the order of a typical quark or lepton mass. Diagonalization 
yields eigenvalues mi ^ — m^j/m/j and 1712 ~ rriR. The matrix that converts from eigenstates 
with positive mass to flavor eigenstates is 



(15) 



i rriD/mji 
-imD/rriR 1 

In this scheme, neutrinos that we know about are much lighter than other leptons because 
of the existence of other very heavy neutrinos. Even without heavy neutrinos, the fact 
that the only dimension-5 neutrino-mass operator in standard-model fields is of Majorana 
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form (jWeinbergl . 



19791 ) leads one to expect that massive neutrinos will be Majorana particles 



In the general iV-fiavor case 

N N 

m = E ^ikPKpk , ^iR = J2 ^I'kPR^k , (/' = / + iV) , (16) 

k=l k=l 

where Pl and Pr are projection operators onto states of well-defined chirality. We as- 
sume here that something like the see-saw with very large m/j is operating so that half the 
eigenstates are very heavy. In that case the light eigenstates mix nearly entirely among 
themselves, and the N x N matrix U responsible, defined to be the "upper-left quarter" of 
U, is nearly unitary: 

N 

UlL ^ E UlkPL(l>k . (17) 
k=l 

Although we have used the see-saw to arrive at Eq. [T71 a mixing matrix U can be defined 
even if the right-handed sector is light or absent. 

The matrix U, which we introduced in Eq. ([3]), nominally has A^^ parameters, A^(A^ — 1)/2 
angles and N{N + l)/2 phases. of the phases are unphysical, so that there are A^(A^ — 
l)/2 independent physical phases. For three active neutrino flavors, the three-phase mixing 
matrix can be written in the form 



U 



iS \ 



-S12C23 - Ci2S23Si3e"^ C12C23 - Si2S23Si3e*'' 



diag{e~, e~ , 1} 



;i8) 



C12C13 S12C13 sise 

S23C13 

iS 

\ ■S12S23 — Ci2C23Si3e — C12S23 — Si2C23Si3e C23C13 J 

where Sij and Cij stand for the sine and cosine of the angles 6ij, 5 is a "Dirac" phase analogous 
to the unremovable phase in the CKM matrix, and the other two phases ai and 02 affect only 
Majorana particles. If neutrinos were Dirac particles these two phases could be absorbed 
into redefinitions of the fields. The M ajorana mass terms in the Lagrangian, however, are 



not invariant under such redefinitions. 



Kobzarev et al 



(119801 ) gives a detailed discussion of 



the number of free parameters in the mixing matrix. 
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III. RATE OF DOUBLE-BETA DECAY 



The neutrino masses and mixing matrix figure prominently in neutrino-mediated j3j3{Qu) 
decay. The rate for that process is 

[^hr = E / \zo.mEe. + - Q,/^^^ , (19) 

spins 

where Zq^, is the amphtude and (5/3/3 is the Q- value of the decay. The amplitude is sec- 
ond order in the weak interaction and depends on the currents in the effective low-energy 
semileptonic Hamiltonian (Ti/?), which we assume for the time being is purely left-handed: 
H/3{x) = G F / V^{e{x)'jf_i{l — 75)z/e(x) Jj^(x)} + h.c, with the charge- changing hadronic 
current. We assume as well that only the particles we know about exist, or that any others 
are too heavy to measurably affect /3/3(Oz/) decay. Then the decay is mediated solely by the 
exchange of three light neutrinos and the amplitude contains a lepton part (a function of x 
and y to be contracted with a similar hadron part and integrated over x and y) of the form 



E e(x)7^(l - 75)f/efc0fc(a;) e(2/)7^(l - i^)Uek(t)k{y) = 

k I 1 

- E ^(a;)7/,(l - i5)Uek(l)k{x) (pl{y)'y^{l + -i^)Ueke''{y) 

k I 1 

where the underbrackets indicate contraction. With our convention 61 = 



(20) 



')k, the contraction 



of (j)k with 0^ turns out to be the usual fermion propagator, so that the lepton part above 
becomes 



47 V (27r)4 



e(a;)7^(l -75) 



.2_^2 ^^y^ + 75)e [y] U^k , 



21) 



k v-'v q^-m 

where q is the 4-momentum transfer. The term with vanishes and the rrik in the denom 
inator can be neglected for light neutrinos, so that the amplitude is proportional to 



("^/3/3) 



k 



nil 



(22) 



The absolute value has been inserted for convenience, since the quantity inside it is squared 
in equation f|T9l) and is complex if CP is violated. Applying the first expression in Eq. fl22|) 
to our one-flavor example, one can see explicitly that a Dirac neutrino, which is equivalent 
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to a degenerate Majorana pair 0i and 02, cannot contribute to /3/3(Oz/) decay because the 
two states would have Uei = 1 and Ue2 = ^, as in Eq. f|T4|) . 

To complete the calculation, one must multiply the lepton part of the amplitude by the 
nuclear matrix element of two time-ordered hadronic currents and integrate over x and y. 
We can write the matrix element of a product of currents between initial and final nuclear 
states i and / as 

n 

where the |n)'s are a complete set of intermediate nuclear states, the -En's are the corre- 
sponding energies, and Ei and Ef are the energies of the initial and final nuclei. When the 
two time coordinates, xq and yo, are ordered and integrated over, and when the exponential 
in Eq. (1231) . is combined with a similar factor from the lepton currents and the g-dependence 
of the neutrino propagator, the following factor in the amplitude results: 



27T5iEf + Eel + E,2 - ^) E 



{fmx)\n){nmm ^ {fmx)\n){n\J^{y)\t) 



(24) 



_q%En + gO + E,2 - Ei) q%En + + E,i - Ei 

We have ignored the neutrino masses, which are small compared to their momenta. The 
quantity q^ = q is the energy of the virtual neutrino, to be integrated over along with the 
virtual momenta, and E'ei, Ee2 are the energies of the outgoing electrons. The energy go 
is typically about an average inverse spacing between nucleons, 100 MeV or so. This value 
is much larger than the excitation energy of states contributing to the decay amplitude, so 
the intermediate-state energies are usually replaced by an estimate E of their average value. 
Studies show that the resulting amplitude is not very sensitive to the choice of the average 
( jPantis and Vergadod . 11990), and that the error it causes in the decay rate is perhaps 15%. 
In this "closure approximation" one replaces the (now unweighted) sum over intermediate 
states by 1, so that Eq. (p^ above becomes 



2'K5{Ef + Eel + Ee2 ' Ei 



(25) 



_q^{E + gO + - Ei) q^{E + gO + - E, 

To go any further, we need an expression for the hadronic current Jl{x). In "the impulse 
approximation" , the hadronic current is obtained from that of free nucleons, and the resulting 
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one-body operator Jl(x) = Y.aOa{x)T^ (where the operator Oa{x) acts on space and spin 
variables of the a*^ nucleon) is used to evaluate the matrix element between states in the 



)Jr(x) because (r. 



+ \2 



0. 



initial and final nuclei. In this approximation Jl{x)J1{y 
The charge-changing hadronic current for a nucleon [i.e. the proton-bra neutron-ket matrix 
element of the current) is 



{p\J^{x)\p') = e^^Mp) [gviq^)!' - 9A{q'Hr - '9M{q^)^q. + gp{q')l,q'^ W) 

(26) 

where q = p' — p, gy = gvi^) = ^, dA = fi'A(O) = 1.26, conservation of the vector current 
tells us that gM^q"^) = gMdviq^) with gM = 5'a/(0) ~ 3.70gv, and the Goldberger-Triemann 
relation, accurate enough for our purposes, says that gp{q^) = 2mpgA{q^) / {q^ + "^^), with 
TUp and m^r the proton and pion masses. The behavior with of the other coefficients can 
be parameterized in a simple way from experimental data: 



9v{q^ 



9v 



[l + qV^VY 



9A{q^ 



9A 



2 \2 



(27) 



with Ky = 0.71 (GeV)^ and = 1.09 (GeV)^. After reducing to nonrelativistic kinematics 
to obtain 0{x), keeping all terms of 0{l/mp) except "recoil" terms that depend on p + p' 
rather than q, and omitting 2"'^ order terms in Ei^2/qo, which takes the electron energies 
out of the denominators of Eq. (!25l) . one can integrate the rate over electron phase space 
making the long-wavelength approximation. Only 0+ 0^ decay is considered so that the 
electrons are predominantly in s states and the effect on them of the charged nucleus as 
they exit may be approximated via a Fermi function. The rate then takes the form given in 
Eq. (ED: 



Oi/ 1-1 
/2\ 



where Qpp = Ei — Ej, Gou{Qi3i3, Z) comes from the phase-space iri t egral, which inc 



Z-dependent Fermi function, and according to 



Rodin et al. 



mm -. 



Simkovic et al. 



(28) 



udes the 



mm . 



M, 



(/I 



2R 



t^9a 



7>j 

qdqJ2 



joiqrab) [hpiq) + hGT{q)^a ■ ^b] + 32{qrab)hT{q) Wj ■ Tab^k ■ Tab " " ^b] 



a,b 



q + E-{Ei + Ef)/2 



(29) 
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Here, the nucleoli coordinates are all operators that, like spin and isospin operators, act 
on nuclear wave functions. The nuclear radius, R, is inserted to ma ke the inatrix element 
dimensionless, with a compensating factor in Gq,^. (As pointed out in ICowelll (120061 ). errors 
have resulted from using different values of R in Mq^, and Gqu.) The internucleon position 



vectors are defined by r, 



ab 



r^l and fa 



rb)/^ab, while jo and j2 are spherical 



Bessel functions, and 



hF{q) ^ -gUr) (30) 

hoAq) = 9a{<1 ) + + ^-r- (31) 

, . ^ _ 9A{q^)9p{q^)q^ 9p{(i^)q^ , 9lAq^)q^ /oo^ 
^^^^^ = + • ^^^^ 

The terms above containing qm are negligible, but those with gp typically reduce the 
matrix element by about 30%. In most calculations, however, even these terms are neglected, 
so that the matrix element takes the approximate form 

Mo. - M^J - 'iMi (33) 
9a 

with 



Mi = {f\Y.H{rab,E)r+T^\{) , and (34) 

a,f) 

Mg^ = {f\J2H{rab,E)aa • a,r+r+|z) . (35) 

a,b 

Here the "neutrino potential" , H is defined as 

— , 2R r°° , sin or 
H(r, E)^ — / dq = — --^ — . (36) 

^ ' Tir Jo ^q + E - {Ei + Ef)/2 ^ ' 

For later reference, we also give an approximate expression for the rate of I3(3{2v) decay, 
which doesn't depend on neutrino mass or charge-conjugation properties, and involves no 
neutrino propagator: 

[ry,]-i = G2.(Q,,, Z)\M§J - ^^Mif (37) 

9a 
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where= G2u{Qi3i3, Z) is another phase-space factor, presented earher in Eq. ([T]), and 



- \ K-(M. + M,)/2 

mGT _ {f\T.a^art\n){n\Y.b^br^\i) /oq\ 
~ V E^-{M, + Mj)/2 ■ ^^^^ 

Nearly all the Fermi strength goes to the isobar analog state in the daughter, so that 
can be neglected. 

We know that there are three light neutrinos with largely left-handed interactions, so it 
makes sense to calculate the /3/3(0z/) rate that those neutrinos and interactions induce. But 
most theorists believe that unobserved particles and interactions exist as well. The most 
popular explanation of small neutrino masses is the see-saw mechanism, which implies the 
existence of heavy neutrinos that couple to left-handed gauge bosons. One simple extension 
of the standard model that gives rise to a see-saw is the "left-right" symmetric model, in 
whic h a heavy right-hand e d weak boson Wpi coe xists alongside the familiar, and lighter, W^. 



Refs. 



Hirsch et al. 



1996 



1 



Prezeau et al. 



(120031 ) give a general analysis of double-beta decay 
in such models. We will not repeat that here, but instead examine the general question of 
whether we can expect physics beyond left-handed weak interactions and light Majorana 
neutrinos to generate double-beta decay at a level that competes with Eq. ( l28l) . 

Right-handed currents can cause /3/3(0z/) through the exchange of both light and heavy 
neutrinos. The coupling of Wr to neutrino mass eigenstates contains a factor Uin (where 
/' labels the right-handed states with definite flavor), while the coupling of the usual 
contains Uu (where I labels the left-handed states), so that the exchange of light neutrinos 
with a right-handed W involved is proportional to 



^ rrikUkUik. (40) 

fc=light 

As we see in our one- flavor example, Eq. ( |T5l) . this quantity is ~ mknici/mR and the am- 
plitude is very suppressed. The largest contribution, not including the one proportional to 
{mpp) derived above, generally comes from the exchange of heavy neutrinos, through two 
WrS. Then there is no suppression from the mixing matrix, the neutrino propagator is 
roughly proportional to l/niR and, crudely speaking. 
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^heavy f / 1 \ , . 

For light neutrino exchange, there are no IVr's in the dominant term, and the propagator 
is roughly proportional to (m^/3)/(g)^, where {q) ~ 100 MeV is a typical virtual-neutrino 
momentum. Then, instead of Eq. ( |4TI) . we have 

flight Gl{mf3i3) 

2(g)2 ^^^^ 

so th at the two amp l itude s will be approxim ately equal when (assuming that Mw^^ ~ 



rriR) (jCirighano et aU . 12004 : Mohapatral . Il999l l 



5 



mn^l'-T^] ' (43) 



which is on the order of 1 TeV for {mpp) y^Am^^^^. Thus, if the heavy mass scale in left- 
right symmetric models is about a TeV or less, it will not be so easy to determine the mass 
scale of the light neutrinos from double beta decay. The same statement is true of many 
other kinds of hypothetical lepton-number-violating models (super symmetry, leptoquarks, 
. . . ) because they usually generate double-beta decay in a similar way, through graphs in 
which heavy particles of some kind play the role of the Wrs and heavy neutrinos. 

Neutrinoless double beta decay in extra-standard models gives rise to new nuclear matrix 
elements. The presence of a single right-handed lepton current causes the g^7p term in 
the propagator of Eq. (I2T]) to contribute to the amplitude, giving rise to derivatives of 
the neutrino potential presented here or forcing one of the electrons into a p state. The 
outgoing p wave leads to a different dependence on the angle between the two emitted 
electrons that could in principle be exploited to distinguish between the action of right- 
handed currents and the neutrino mass in light neutrino exchange. But the short-range 
exchange of a heavy particle will not always manifest something like the g''7p term, and 
often the only way to distinguish such a process from neutrino-mass-induced decay is to 
exploit the different nuclear matrix elements that enter. Provided the matrix elements can 
be accurately calculated, analysis of measured lifetimes in several isotopes or to several 
states in the same istotope can tell you whether long or short range is responsible (a little 
more on this below). Of course, as already mentioned, the accuracy with which nuclear 
matrix elements can be calculated is a big issue, and we discuss it later. A more detailed 
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treatment of the matrix e 
found in Refs 



mm; 



Doi et al. 



Tomodal ( 



199l|). 



emen t s governing the various kind s of double-beta decay can be 



(119851); 



Haxton and Stephenson Jr. 



fll984h : 



Simkovic and Faessler 



The imphcations of some popular extra-standard models for (3(3{Qv) are discussed below. 
We close this section with two general points. First, when lepton number is spontaneously 
broken, as it is in most models that result in a see-saw mass matrix, there must exist one or 
more zero-mass bosons, called Majo rons, that could be emit t ed along with the two electrons 



in double beta dec ay (/9/3(0z/, x)) (jChikashige et al. 



1981 



Gelmini and Roncadelli 



1981 



Georgi et ali Il98ll ). Apparently, however, it is difficult for such a process to have a very 



large amplitude. Second, even if some exotic lepton-number-violating physics exists and light 



neutrino exchange is not responsible for the decay, t he occurrence of /3/3(O z /) sti 



1 implies that 



neutrinos are Majorana particles with nonzero mass (iSchechter and Valid . 119821 ). The reason 
is that any diagram contributing to the decay can be inserted into a neutrino propagator, 
with outgoing electron fines closed appropriately as in Fig. [3l If /?/5(0i/) decay is observed, 
we will know for certain that neutrinos are their own antiparticles, even if the possibility of 
exotic physics or uncertainty in the nuclear matrix elements prevents an accurate extraction 
of the neutrino mass scale from the observation. 



IV. DOUBLE-BETA DECAY AND NEW PHYSICS 



Over the past few decades, especially the last, much has been learned about the neutrino 
mixing angles and mass eigenvalues. Table [T] summarizes our knowledge of these neutrino 
parameters. These results have increased the importance of /3/3(0i/) experiments; in the first 
subsection below, we explain why. The other subsections discuss other physics that might 
be revealed by (3[3{Qv). 



A. Neutrino Mass 



If neutrinos are Majorana particles they will mediate f3(3{0i>) at a rate proportional to the 
square of {mj^p), Eq. 0221) . The known values of the mixing- matrix elements in Eq. fll8p allow 
us to predict the rate of /?/3(0z/) under several scenarios for the neutrino's mass spectrum. 
If we ignore the LSND result (see Sec. flIV.C|l ) the oscillation data are consistent with only 
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3 such masses, but their spectrum can still take 4 possible forms: 

1. Normal Hierarchy-Dirac The two masses with the smaller splitting indicated by Sm'^^i 
are smaller than the 3^'^ mass. The neutrinos are Dirac. 

2. Inverted Hierarchy-Dirac The two masses with the smaller splitting indicated by Sm'^^i 
are larger than the 3'"'^ mass. The neutrinos are Dirac. 

3. Normal Hierarchy-Majorana The two masses with the smaller splitting indicated by 
Sm^^Y are smaller than the 3^"^ mass. The neutrinos are Majorana. 

4. Inverted Hierarchy-Majorana The two masses with the smaller splitting indicated by 
Sm^^i are larger than the 3'^'^ mass. The neutrinos are Majorana. 

In addition, since the absolute mass scale is unknown, it is possible that the differences be- 
tween the 3 mass eigenvalues are small compared to the masses themselves. In this arrange- 
ment the neutrinos are referred to as quasi-degenerate or sometimes simply as degenerate. 
In the normal hierarchy, the state corresponding to the largest mass contributes with a small 
mixing angle. Hence if the mass of the lightest state, mughtest is small, is also small. 

By contrast, in the inverted hierarchy, the heavier neutrinos are large contributors to {mpj^). 

Modern (3(3 research is exciting in part because if nature has selected possibility 4 we 
should be able to see (3(3{Qv) with the next generation of experiments. By the same token, 
null experiments would rule it out, and could restrict the parameter space available for 
possibility 3. And if some very sensitive experiment ever saw a very small {nipp), it would 
demonstrate that possibility 3 is nature's choice. Actually certifying that possibility 4 is 
the choice is a trickier matter, though. To see this, we show in Fig. H] the dependence of 
(mpis) on miightest- To make this plot we have used the best fit values of the mixing angles 
(^soi = 33.9°, ^atm = 45° and ^13 = 0°) and the Am^'s (Table ^. The figure shows, as 
expected, that {m/sp) is larger in the inverted hierarchy than in the normal one. (The plot 
shows regions rather than lines due to the unknown Majorana phases.) But because there is 
no way of measuring mughtest, PPi^^) will not be able to distinguish the inverted hierarchy 
from the quasidegenerate arrangements. (Although in principle it is possible to directly 
measure mughtest in a beta decay experiment with ideal energy resolution, in practice it is 
not feasible.) A large {frtpj^) will still not tell us for sure which eigenstate is the lightest. And 
we also won't know for sure that other TeV-scale physics isn't responsible for the decay. 
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Other measurements can help, however. Unhke /3/3(0z/), the rate of which reflects the 
coherent exchange of virtual neutrinos, beta decay involves the emission of real neutrinos, 
whose mass can alter the beta particle spectrum. The corresponding effective beta-decay 
mass {mji) reflects the incoherent sum of the mass terms: 



(^/3) = E ^3 \Uej\ = ml\Uei\^ + ml\Ue2? + m2|?7e3p . (44) 
-j 

In a beta-decay experiment this quantity would approximate the difference between the 
endpoint of the electron spectrum and the Q-value. The approximation is valid as long as 
the energy resolution is too poor to separate individual endpoints due to each of the m^. 
For all presently planned beta decay experiments, that is the case. 

Equations (12^ and (jH]) depend differently on the mixing angles, phases and the mass 
eigenvalues. If beta-decay experiments find {mp), PP^Ou) measures {mpj^), and Mq^ ever get 
accurate enough, they could help constrain the Majorana phases discussed below. If (m^) 
yielded a result that was inconsistent with the two known Am^'s and a measured {mpjs), it 
could demonstrate that new physics, either in the form of new particles exchanged in /5/5(Oz/) 
decay or sterile neutrinos that mix with the 3 active neutrinos, is at play. 

We should no t e tha t the neutrino contribution to the mass density (fij,) of the uni- 



verse ( iHannestedl . l2003l ) constrains a third combination of the neutrino masses: 



n^h^ = , (45) 

92.5 eV ' ^ ^ 



where 



T, = rrij = mi + 1712 + . (46) 

3 

Since no experiment measures the mass eigenstates directly, effective neutrino mass measure- 
ments, coupled with the oscillation measurements, are all required to determine a complete 
set of the parameters describing neutrinos. 



B. The Majorana Phases 

The elements of the mixing matrix may contain many phases, but only a small number 
have any physical consequence. In a 3 x 3 mixing matrix describing Dirac neutrinos, all 
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but one phase (the so-called Dirac phase S) can be absorbed into a redefinition of the 
neutrino and anti-neutrino fields. If neutrinos are Majorana, the phases of each neutrino 
and antineutrino pair are correlated and fewer phases can be eliminated in this way. For the 
3x3 case, two additional phases (the Majorana phases ai, ai of Eq. flTSl) ) are allowed. Any 
of these phases can result in CP violation if its value differs from a multiple of vr. In Fig. HI 
the borders of the shaded regions correspond to phases that are multiples of vr, resulting in 
CP conservation. The shaded regions between the borders correspond to CP violation. By 
comparing measurements of S or (m^) to (mpf^), one might be able to demonstrate that CP 
is violated in the lepton sector. 

The Dirac phase can lead to observable effects in oscillation experiments, whereas the 
Majorana phases that appear in Eq. fl2^ have no effect in those experiments. To see this, 
note that the transition amplitude for oscillations from a neutrino of fiavor / to fiavor j is 



Ziui 



m, L 



(47) 



where L and E are the distance traveled and energy respectively of the u. The Majorana 
phases in the diagonal matrix of Eq. f|T8|) cancel with their complex conjugates in the product 
UU*. 



Pascoli and Petcovl . 



2003 



Sugiyama 



2003f ). lEUiott and Engell ( 120041 ) provided a graphi- 



i 


■ pre- 




2003; 



cal example of how the measurements from oscillations, /3/?(0z/), beta decay, and cosmology 
might be combined to learn about possible values for the phases of Eq. |22l But any attempt 
would have to contend with the problem that there are two such phases but only one known 
experimental measurement that is sensitive to them. With only /5/5(Oi/) as a probe, no 



unambiguous measurement o 



violation in the lepton sector (iBuchmuUer et al. 



both g's is possible. A 



2002 



though Leptogenesis depends on CP 



Fukugita and Yanagida 



19861), it will 



not be easy to use cosmological measurements to help quantify the a's, because the relation 
between CP violation in the decay of heavy Majorana neutrinos in the early Universe and 
the phases in the light-neutrino mixing matrix is model dependent. 

If 6i3 is zero, however, only one of the a's will contribute to (jripp). This might allow 
that a to be determined from measurements of /9/3(0z/). There is still a problem, though: 
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extracting information about a will require the ca l culati on of M^y with an accuracy that is 
currently out of reach (IBarger et alX . l2002l : lEUiottl . l2006l ) . 



C. Sterile Neutrinos 



The LSND (lAguilar et a/.l . l200ll ) neutrino-oscillation result indicates a Am^ scale of ~1 



eV^. A value this large cannot be incorporated into a 3-neutrino mixing scheme along 
with the atmospheric and solar scale Am^ values. This inconsistency is sometimes called 
the Three I\vp? Problem: with 3 neutrino masses there are only 2 mass differences. One 
approach to solving this problem is to add some number of light sterile neutrinos to the 
three active neutrinos known to exist. 

Lig ht sterile neutrin o s might seem unlikely, but can be produced by the see-s aw mecha- 



nism ( 1 Cell- Mann et al. 



1979 



Mohapatra and Senianovid. 



metry makes t 



l e Ma i orana mass matrix s i ngula r ( 



Goldman et al 



2000 



Stephenson Jr et all l2005l ). 



1980 



Yanagidal. 



Chikira et al. 



2000 



19791) if a sym 



1998 



Chun et al. 
Goswami and RodejohanJ ( 20061 ) have 



reviewed the consequences of light steriles of yarious kinds for j3f3. In particular they discuss 



a model, originally suggested by 



Sorel et al. 



( 120041 ) that adds two sterile neutrinos. One 



solution reproduces the LSND data while still matching the null results from other short- 
baseline-accelerator neutrino oscillation experiments and does not violate constraints from 
cosmology. This solution (Am^^ = 0.46 eV^, Am^^ = 0.89 eV^, Uei = 0.090, Ue5 = 0.125) 
would provide a maximu m additional contribution to (m/^/^) . " ^ ^0 m eV if mi is small. The 



MiniBooNE experiment (lAguilar-Arevalo et al 



2007 



Bazarkd . 



2000l ) has found no evidence 



for the 2-fiavor oscillation model as an explanation for the LSND result, although the exper- 
iment sees a low-energy excess that is unexplained as of this writting. Whether models with 
3 active neutrinos and 2 or more sterile neutrin os can explain both sets of data is still being 
investigated {e.g. iMaltoni and Schweta (120071 )). Sterile neutrinos, however have the virtue 
of saving the hea vy element nucleos ynthesis process in the neutrino-driven wind behind a 
supernova shock (IFetter et ali 120031 ) . and are not ruled out in general by any experiments. 
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D. Alternative Mechanisms 



If /3/3(0z/) is observed, it will, as discussed in Sec. IIIIl demonstrat e that the neutrino 



mass matrix has a nonzero Majorana component (iSchechter and Valid . 119821 ). But, as we 
also discussed earlier, this does not imply that the decay rate is proportional to the square 
of {nif^ij), as it would be if the light-neutrino mass were the driving mechanism. After any 
observation of /3/5(Oz/), thefore, it would become important to find an experimental signature 
that isolates the underlying physics. The decay rate for PP^Ou) depends on a Lepton Number 
Violating Parameter (LNVP) — (^/j/j) in the case of light neutrino exchange produced by 
left-handed currents — that will be different for each possible mechanism. In addition, the 
nuclear matrix element may also depend on the exchange mechanism. 

A number of suggestions for potential signatures of the underlying physics have been 
published. As noted above, if the weak interaction includes right handed currents, the 
left-handed virtu al neutrino could couple to the absorbing neutron without an helicity flip. 
Doi et all ( 119851 ) suggested the use of the kinemat ic distributions to discern right-handed 
currents from light-neutrino- mass contributions and lAli et all (120061 ) discussed the use of the 



angular distrib ution to distinguish the left-right symmetric model from other possibilities. 



Tomoda 



2000l ) proposed examining the ratio of rates to the 2+ excited state to that to the 



ground state in the same isotope as a signature of right-handed currents. 

Many alternative exchange mech anisms invo lv e heav y exchange par ticles ( ^ 



TeV), such as heavy neutrinos [ 



1996b), supersymmetric particles ( 



Hirsch et al. 



Hirsch et al. 



ears (iKlapdor-Kleingrothaus and Sarkai 



1996af). lep t oquar ks (IHirsch et al. 



1996c 



Vergadod . |2002| ) and scalar bilin- 



20031 ). As noted alread y, electron-kin e matic s 



will not in general d epend on the heavy particle exchanged (jPrezeau et all |2003| ). 



Simkovic and Faesslerl (120021 ) showed that the relative rates of the decay to the first ex- 



it-neutrino exc hange, the 



Prezeaul (120061) estimated 



cited 0+ state and the ground state might distinguish among ligl 
heavy-neutrino exchange, and supersymmetric-particle exchange, 
the relative contributions to P^'^ from light neutrinos and heavy particle exchange. 

Some models lead to complic ated expressions for decay rates. Singlet neutrinos in extra 



dimensions can mediate /3/?(0i/) (IBhattacharyya et al. 



20031 ) , but the Kaluza-Klein spectrum 



of neutrino masses spans values from much less than to much greater than the nuclear 
Fermi momentum. One cannot, therefore, factor the decay rate into a product of nuclear 
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and particle factors. Futhermore, the decay rate depends on unknown parameters such as 
the brane-shift parameter and the extra-dimension radius. Finally, mass-varying neutrinos 
(MaVaNs) might lead to decay rates that depend on the density of the matter in which the 



pr ocess occurrs (jKaplan et al. 



200J). 



Cirigliano et al\ (120041 ) recognized the potential of using /i — e and |U — 67 in combi- 



nation with j3(3{Qv) to decide whether light-neutrino exchange is the dominant mechanism. 
Certain supersymmetric particles and heavy Majorana neutrinos could be produced at the 
Large Hadron Collider. Seeing these particles, or not seeing them, could help us determine 
what is responsible for /3/3(0i/). 



V. CALCULATING NUCLEAR MATRIX ELEMENTS 



A. Nuclear-Structure Theory 

The better we can calculate the nuclear matrix elements that govern double-beta decay, 
the more accurately we can extract an effective neutrino mass, or a heavy-particle mass. 
Some scientists have tried to minimize or play down the uncertainty in the matrix element 
in an attempt to strengthen claims regarding the precision with which {mpp) can be deter- 
mined. In contrast, physicists in other fields have occasionally exaggerated the uncertainty 
to downplay the impact /5/?(0z/) would have regarding conclusions about neutrino mass. Here 
we try to assess this complex issue as objectively as possible. To begin, we must review the 
accomplishments, difficulties, and main lines of inquiry in the theory of nuclear-structure, 
which is currently in an exciting period. 

There are a large number of nuclei and many phenomena to observe in each one. Nuclear- 
structure theory has tended to concentrate on properties that are easiest to measure: spectra, 
electromagnetic transition rates, and cross sections for the scattering of electrons or hadrons. 
Traditional nuclear-structure theory divides the chart of nuclides into different regions — 
some in which the nuclei are deformed, some in which they are spherical, others in which 
the behavior is more complicated — and adopts different models for each kind of nucleus. 
Increased computing power is allowing us to gradually abandon these models, which usually 
consider the dynamics of only one or a few nucleons explicitly, lumping the rest into a 
collective core with only a few degrees of freedom. Calculations now are increasingly "a6 
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initio" , attempting to solve the many-body Schrodinger equation directly from empirical two- 
and three-nucleon interactions, which are obtained from fits to A^A^ phase shifts and energy 
levels in the triton and ^He. Monte Carlo techniq ues now yield accura t e wav e functions with 
such interactions in nuclei as heavy as ^^C (see iPieper and Wiringal (120011 ) and references 
therein), and controlled approximation sch emes that treat all th e nucleons with the same 



interaction s (e.g. the no- c ore s. 



approach (IKowalski et al. 



lell model (^ 



2004 



Wloch et al 



avrati 



et al. 



2OO0I ) and the coupled-clusters 



2OO5I )) currently work well for bound-state 



energies in nuclei up to about ^^O, with the closed-shell nucleus ^°Ca in range. 

Unfortunately, these calculations aren't likely to reach atomic numbers relevant io (5 j3 
soon, and even if they did, would focus first on spectra and other observables generated 
by one-body operators before considering the difficult two-body decay operator discussed in 
Section IIIII The other main line of structure theory, however, is geared towards systematic 
(though less accurate) predictions in arbitrarily heavy nuclei. Its framework is mean-field 
theory (treating all the nucleons) and extensions such as the random phase approximation 
(RPA), restoration of symmetries artific ially broken by the mean field, and the generator- 



coordinate method ( IBender et all l2003l ). These techniques, which can involve modifying 
the equations in ways that improve phenomenology but at first sight appear dubious, are 
closely related to density-functional theory, which leads naturally to similar mean-field-like 
equations and has been spectacularly successful in electronic systems. But as in lighter 
nuclei, the initial focus is on nuclear spectra, with some attention to simple transitions. 
Though other observables can be calculated, it can be difficult to judge the accuracy of the 
results, because systematic data with which to test calculations is not available. 

These two lines of research — - ab initio and mean-field calculations — are both developing 
rapidly and naturally absorb most of the energy of the nuclear-structure-theory community. 
Relatively few researchers focus on double-beta decay, which involves complicated operators 
in structurally complex nuclei such as ^^Ge, and for good reason: there are so many more 
general problems that have frustrated theory for years but are now becoming solvable. As 
a result, the models and techniques that have been applied to double-beta decay are more 
limited in scope and power than they might be. They fall primarily into two categories, the 
Quasiparticle Random Phase Approximation (QRPA) and the Shell Model. We discuss each 
below and try to assess the accuracy of the corresponding results before looking at prospects 
for improved calculations. 
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B. QRPA 



The neutron-proton QRPA, first used in iHableib and SorensenI (119671 ). is an extension 
of the random-phase approximation, a technique primarily for calculating the energies and 
excitation strengths of collective states. In calculations of double-beta decay, a number of 
approximations are made even before the QRPA is applied. Most nuclear-structure calcula- 
tions add to the Hamiltonian, a one-body potential/mean-field, which is diagonalized first to 
provide a single-particle basis in which to treat the two-body (and in principal three-body) 
Hamiltonian that comes from A^A^ and NNN data. In a complete calculation, the effects 
of the one-body potential are subtracted off later. In most nuclear-structure calculations, 
however, the potential is retained and the single-particle basis truncated. In double beta- 
decay QRPA calculations, the "model space" typically includes states within 10 or 20 MeV 
of the Fermi surface. The energies of these states (or equivalently, the one-body potential) 
are adjusted to reproduce properties of states with a single nucleon outside a closed shell. 

The two-body interaction is also not treated exactly, again even before application of 
the QRPA. It is possible in principle to construct an effective Hamiltonian and effective 
transition operators that, for low- lying states, have the same eigenvalues and matrix elements 
in the truncated single-particle space as does the full Hamiltonian in the full Hilbert space. 
In appli cations of the QRPA, the effective interaction is usually approximated by a G- 
matrix (IHjorth- Jensen et all Il995l ) , which takes into account the exclusion of very high- 
lying levels but not those just above the model space or those below, which are assumed 
to contain an inert core. Then certain parts of the interaction and operator are scaled 
to reproduce data: all 0"^ matrix elements are sometimes multiplied by about 1.1 so that, 
when the interaction is treated in the BCS approximation (the starting point of the QRPA) 
they reproduce experimental odd-A — even-A mass differences. The parameter that is 
squared in the double-beta operator is often taken to be 1.0 rather than 1.26, and short- 
range correlations responsible for converting the bare Hamiltonian into a G-matrix must be 
inserted explicitly (although approximately) into the operator's two-body matrix elements. 
Other more dramatic adjustments take place in the application of the QRPA itself. 

To understand the logic behind the neutron-proton QRPA, we first remove the "Q", 
addressing the simpler charge- changing RPA. The two-body operator that governs (in 
the closure approximation) the decay, in second-quantized form, can be written O = 
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J2ijkiOijkiPiPjiT-kn-h where p' creates a proton, n destroys a neutron, and the coefficients 
Oijki are two-body matrix elements. One can rewrite this operator as a sum of products of 
one-body operators O = J2ijkiOijki{plnk){plni). In the charge-changing RPA the one-body 
p'^n and n'^p operators are replaced by operators with the same quantum numbers but boson 
commutation rules, both in the decay operator and the Hamiltonian. After the "bosonized" 
Hamiltonian is diagonalized to obtain a representation of the intermediate-nucleus states 
as a set of "one-phonon" states, the bosonized one-body operators contained in the decay 
operator then connect those states with the mother and daughter ground states, which are 
"phonon vacua". (Actually there are two sets of intermediate states, one corresponding 
to excitations of the mother vacuum and one to excitations of the daughter. The sets are 
not identical, a fact that has so far been taken into account only approximately.) This 
"quasi-boson" approximation is expected to work best for strongly excited intermediate 
states, where excitations with different single-particle quantum numbers contribute coher- 
ently. When the intermediate states are summed over, the result is the ground-state to 
ground-state double-beta decay matrix element. 

The QRPA differs from the above only in that pairing between like particles is taken 
into account through BCS theory, so that the mother and daughter nuclei, are now "quasi- 
particle phonon" vacua. They include correlations by breaking number conservation and 
have the correct number of neutrons and protons only on average. The smearing of the 
Fermi surface induced by pairing has a large effect on the matrix elements, which actually 
vanish for /9/9(2z/) in the absence of pairing because the transition to the daughter ground 
state is Pauli blocked. 

Two-quasiparticle states, the creation operators for which are now the bosonized objects, 
contain some components that are particle(proton)-hole(neutron) excitations and others 
that are particle-particle or hole-hole excitations (based on nuclei with two more or less 
particles than the one we're actually interested in). The interaction has matrix elements 
in both channels, and these matrix elements are usually treated as independent because a 
large change in, e.g., a single particle-particle matrix element translates into small changes 
in all the particle-hole matrix elements. The particle-hole matrix elements are generally 
left alone because their primary effect is on the giant Gamow- Teller resonance, which is 
collective and reproduced pretty well by the QRPA. The particle-particle matrix elements, by 
contrast, affect low-lying states with much less strength, and are multiplied by a parameter 
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traditionally called gpp that takes values between 0.8 and 1.2. This parameter turns out 
to have a larg e effect on double- b eta decay, particular l y the 2i/ variety and particularly 



in the QRPA (lEngel et al. 



1988 



Vogel and Zirnbauerl . Il986l ). It is needed because the 



QRPA simplifies the particle-particle correlations in a way that must be compensated. The 
parameter is usually adjusted to reproduce related observables such as or /5/?(2z/) rates, 
which also depend on it strongly. But although adjusting the parameter increases (3^ rates, 
for example, to better agree with experiment, it appears to concentrate the strength of the 

operator at too low an energy. And if the parameter is increased a bit past its best value, 
an unphysical deuteron condensate becomes lower in energy than the BCS-like quasiboson 
vacuum, causing the QRPA to break down completely. Some of this problem is due to the 
artificial separation of particle-particle correlations into "static" correlations associated with 
isovector pairing that are in the BCS quasiparticle vacuum, and "dynamic" correlations, 
associated with isoscalar pairing {i.e. neutron-proton particle-particle interactions in the 
L = channel) that appear in QRPA through the parameter gpp as a correction to the 
quasiparticle vacuum. A satisfactory treatment of both kinds of pairing on the same footing, 
something we are still awaiting in QRPA-like methods, would help. 

In the meantime, much of the QRPA work has revolved around two questions: 1) How, 
besides the better treatment of pairing just mentioned, can the artificial sensitivity of the 
QRPA to Qpp be reduced, and 2) which observables should fix gpp (beta decay, (3(3{2v) . . .) 
and how accurate will the resulting (3(3{Qv) matrix elements be? The first issue has been 
the subject of an inipressively large number of papers. One a pproach is "Second" QRPA 



(jRaduta et al. 



1991 



Stoica and Klapdor-Kleingrothaud . l200ll ). in which the fact that the 



quasiboson approximation is merely the first term in an exact expansion is exploited by 
calculating the eff ects o f the next term. Another very popular approach is Renormalized 
QRPA (RQRPA) ( Hara . 1964 : Rowe, 1968 : Toivanen and Suhonen . 19951 ). The quasiboson 



approximation is equivalent to replacing commutators by their ground-state expectation val- 
ues, with the ground states treated in the BCS approximation. The RQRPA uses the QRPA 
ground states instead of BCS states to evaluate the commutators. Because the commutators 
in turn help fix the ground states, the two are evaluated self-consistently. A variant of this 
approach is the "Full RQRPA", in which the effects of isovector np particle-particle inter- 
action, artificially strengthened to account implicitly for isoscalar pairing that underlies the 
sensitivity to gpp, are included in the BCS calculation that defines the quasiparticles as well as 
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in the subsequent QRPA calculation (jSchechter et al. 



1996 



Simkovic et al. 



tor np pairing was f irst in troduced in this way in the unrenormalized QRPA (jCheoun et al. 



993 



Pantis et al 



997h. flsovec 



19961 ). Another extension is the Self- Consistent RQRPA (SCQRPA) 



(iBobyk et all l200ll ). in which the occupation numbers in the RQRPA ground state are fed 



back into the BCS approximation and the BCS-RQRPA sequence iterated until the BCS 
and RQRPA ground states are consistent. All these methods reduce the dependence of the 
matrix elements on the strength of the neutron-proton pairing interaction. 

On the other hand, all these calculations also leave out correlations, and the profusion 
of related methods, none clearly superior to the rest, is confusing. Some people have even 
begun to treat the calculations as a statistical sample fr om which to extract an error bar 



on the calculated matrix elements (IBahcall et al. 



2004bt), reachin g the conclusion that the 
matrix elements are only known to within an order of magnitude. Elliott and Engell (12004 ) 
tabulate results for ''^Ge and discusses why a statistical analysis is not a good idea: many 
of the calculations are explicitly prelimi nary or make no attempt to reproduce related data. 



Here, we merely note that recent papers 



Rodin et al 



( 12003 



20061 ). addressing issue 2) above. 



argue that if the strength of the pairing interaction is adjusted to fit experimental mass 
differences and (especially) if Qpp is then adjusted to reproduce measured /3/3(2z/) rates, 
then almost all QRPA-like calculations, independent of their choice of number of single- 
particle levels, interaction, quenching of qa, etc., give the same results to within about 30%. 
If not, then they have either not included the induced-pseudoscalar terms in the nucleon 
current or neglected to add short-range correlations to the wave functions. Though we have 
no way of knowing for sure that any of the calculations give results that are right, it is 
comforting to think that they agree with one another. Not everyone is reassured, however, 
for two reasons. First, the claim that agreement with measured is more important than 
agreement with other observables, such as singl e-beta decay rat es from the lowest lying states 
in the intermediate nucleus, can be disputed (jSuhonenl . l2005l ). When single-beta decay is 
used to adjust parameters, the resulting Po^ are different because the QRPA is not able 
to reproduce both the beta decay and the /3/3(2z/) at the same time. Second, the size of 
matrix-element quenching by short-range correlations, which affect (3(3{Qv) but not (3(3{2v) 
(for which the operator has no radial dependence) , is under debate. Most authors have 



used th e phenomenological correlation function of 



papers (IKortelainen et al. 



2007 



Miller and Spencerl (119761 ). but recent 



Kortelainen and Suhonen 



20071 ) argues that more realistic 
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procedures, e.g. the Unitary Operator Correlation Method (see, e.g., iRoth et al\ (l2005l )). 
produce substantially less quenching. Until these issues are resolved — and the short-range- 
correlations issue clearly can be with more theoretical work — it is hard to say that the 
QRPA gives unique predictions for the matrix elements (and therefore to the value of (mp/s) 
extracted from an experiment) to better than a factor of 2 or so. 

There are also systematic deficiencies in the way the QRPA has been applied, in addition 
to the obvious simplification inherent in the method itself, that may cause errors. Almost all 
calculations so far have assumed that the ground states are spherical, when we know in many 
cases, e.g. ''^G e, that they are not. Some preliminary work on 00(^ 2i>) with a "deformed 
QRPA" exists ( Alvarez- Rodriguez et al . I2OO6 : Isimkovic et aL ,2003, ). but /5/?(0i/) remains 



uninvestigated. With today's computers, however, there is no reason not to use a deformed 
QRPA. The inclusion of the continuum is anoth e r feat ure that has been implemented just 
once, in a preliminary way ( Rodin and Faessler . 20061 ). Here too, the state of the art in 
nuclear structure is more advanced than what is being applied to double-beta decay. In 
mean-field-based nuclear structure studies, continuum versions of the QRPA that treat all 
the nucleons (so that there is no inert core) have been applied many times, a nd versions of 
these treatments that include deformation are beginning or about to appear ( Hagino et all 
20041 ). A Skyrme-HFB-continuum QRPA calculation of /5/3(Oz/) which is clearly possible 
even with deformation included, would be worthwhile. Calculations of single-beta decay 



rates in that f ormulation, a. 



time (see, e.g., 



thou gh not yet with deformation, have been around for some 



Engel et al\ ( 1l999h ). Perhaps one reason this more complete calculation has 



not been done for double-beta decay is that it would still treat isovector and isoscalar pairing 
differently, and as a result would still show sensitivity to gpp. What is really needed to lessen 
reliance on parameter tuning is a more unified treatment of particle-particle correlations. 
The shell model, although it has its own set of difficulties, does provide that. 

C. Shell Model 



The nuclear shel 



model was recently reviewed by expert practitioners in this journal 
( ICaurier et a/.l . l2005l ). Compared with the QRPA it has only one disadvantage: the (much) 
smaller number of single-particle states that can be included. But because the number 
is smaller, one can include correlations of arbitrary complexity within the single-particle 
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space. In addition, it's possible to calculate essentially any spectroscopic observable. In the 
QRPA, by contrast, only energies of or transitions from ground states can really be reliably 
calculated because the correlations necessary to describe configurations more complicated 
than a two-quasiparticle excitation of the ground state are not included. 

Because of its sophisticated treatment of correlations, a good shell-model calculation is 
harder to carry out than a QRPA calculation. First, as in the QRPA, the energies of single- 
particle states, which usually are assumed to have harmonic-oscillator wave functions, must 
be fit to spectra in closed-shell + 1 nuclei and an effective interaction appropriate to the 
model space must be constructed. Typically something like a G-matrix is a starting point 
but with many extra terms subsequently added and fit to a wide range of energies and 
transition rates (painstaking work) in nuclei near the beginning or end of the shell, where 
Hamiltonian matrices are relatively small. In the nuclei of interest, the Hamiltonain matrices 
are much larger — nowadays up to 10^ x 10^ or even larger — and the Lanczos method is 
often used to diagonalize them. The transition operator is modified by inserting the effects 
of short-range correlations, again as in the QRPA. The axial vector coupling constant is 
sometimes quenched, but not always. 

The first modern shell-model calculations of jS/S decay date from 



Haxton and Stephenson Jr.l (119841 ) and references therein. Only a few truly large- 
scale shell model calculations have been performed. The heavy deformed jS/S nuclei, ^^^f/, 
and ^^°Nd, for example, require bases that are too large to expect real accuracy from 
the shell mo del. ^. Realistic work has t hus been restricted to ^^Ca, ^^Ge, ^^ Se, and ^^^Xe 



( ICaurier et al. 



1996 



Horoi et al. 



2OO7I ) (and now, though still unpublished (iPovesl . l2006l ) 



in Cd, Te, and Te as well). Less com prehensive calculations have been carried out 



in several other nuclei (ISuhonen et all 119971 ) . 

The largest-scale shell-model calculations tend to produce matrix elements that are a bit 
smaller than t h ose fr om QRPA, by factors of up to 2 or 3. In addition recent unpublished 
results (jPoveg, 120061) appear to differ by 50% or more from earlier published results by 
the same authors (ICaurier et all Il996l ). who are essentially the only people working on the 
problem in a comprehensive way. Some of the uncertainties, e.g., the quenching of qa and 
the effects of short-range correlations, are the same ones that afflict the QRPA; others are 



^ Psuedo-SU(3)-based truncations have been used for these nuclei {e.g., in Hirsch et al 1996). 
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connected with truncation and determination of the effective Hamitonian. 

The shell model calculations can be improved in the short term, even without go- 
ing to larger spaces ( which is in fact possible via, e.g. the factorization method of 



Papenbrock and DeanI ((200^). Perhaps the most important step that could be taken now is 



a more systematic treatment of the effective decay operator appropriate for the shell-model 
space. Perturbative corrections to the operator that account for the finite model-space size 
can be calculated, alt hough perturbation t heory in the residual two-body interaction does 
not always converge (iHjorth- Jensen et all Il995l ). Its use would nonetheless give us a bet- 
ter idea about the size of corrections to the naive renormalization (quenching of qa and 
artificial short-range correlations) that is the current state of the art. Ex ploratory calcula- 



tions in a solvable model indicate that the corrections may be significant (lEngel and Vogell . 



20041). Though perturb ative corrections to the Gamow- Teller operator have been evaluated 



(ISiiskonen et al. 



200 ll ). nobody has looked at two-body operators beyond the Hamiltonian. 



D. Prospects for the future 

Most good calculations give the same result for a given matrix element to within a factor 
of 2 or 3, usually less. Assuming that there is no many-body physics that affects the result 
but is systematically absent from all calculations, we can take that range to approximately 
represent the level of uncertainty in our knowledge of the matrix element. What are the 
prospects for reducing that uncertainty? 

In the short term we can make progress by being more careful — by quantifying uncer- 
tainty in the weak nucleon current, form factors, short-range correlations, quenching, etc.. 
The size of short-range-correlation effects is currently in dispute, but structure theorists 
now know enough about such effects to resolve the debate. In the medium term, the best 
hope is a better shell-model calculation. Here what we really need is progress in construct- 
ing effective operators, both the Hamiltonian and the decay operator. Recent years have 
seen dramatic progress in our abil ity to compensate for high-momentum physics that is cut 



out (see, e.g., 



Bogner et al 



Dut reliably correcting for low energy excitations such 



as core polarization is a longstanding problem. Partial summation of diagrams, a tool of 
traditional effective-interaction theory, is helpful but apparently not foolproof. 

In the long term these issues will be solved. As already mentioned, the coupled-cluster 
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approximation, an expansion with controlled behavior, is being applied in nuclei as heavy 
as ^°Ca. With enough work on three- and higher-body forces, on center-of-mass motion, 
and on higher-order clusters, we should be able to handle ™Ge. The time it will take is 
certainly not short, but may be less than the time it will take for expermimentalists to 
see neutrinoless double-beta decay, even if neutrinos are indeed Majorana particles and the 
inverted hierarchy is realized. And the pace of theoretical work will increase dramatically 
if the decay is seen. Observations in more than one isotope will only make things better. 
Our opinion is that the uncertainty in the nuclear matrix elements in no way reduces the 
attractiveness of double-beta decay experiments. Given enough motivation, theorists are 
capable of more than current work seems to imply. 



VI. EXPERIMENTAL ASPECTS 



A. Background and Experimental Design 



Double beta decay experiments are searching for a rare peak (see Fig. [5]) upon a continuum 
of background. Observing this small peak and demonstrating that it is truly /3/3(0i/) is 



a challenging experimental design task. 



exper iment have been discussed previously ( lElliott 



he characteristics that make an ideal PPjOiy 



2003 



Elliott and Vogei l2002l : IZdesenkol . 



2OO2I ). Although no detector design has been able to incorporate all desired characteristics, 
each includes many of them. (Section IVII.CI describes the various experiments.) Here we 
list the desirable features: 



The detector mass should initially be large enough to cover the degenerate mass region 
(100-200 kg of isotope) and be scalable to reach the inverted-hierarchy scale region 
1 ton of isotope). 

The /3/5(Oz/) source must be extremely low in radioactive contamination. 

The proposal must be based on a demonstrated technology for the detection of (3 (3. 

A small detector volume minimizes internal backgrounds, which scale with the detector 
volume. It also minimizes external backgrounds by minimizing the shield volume for 
a given stopping power. A small volume is easiest with an apparatus whose source is 
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also the detector. Alternatively, a very large source may have some advantage due to 
self shielding of a fiducial volume. 

Though expensive, the enrichment process usually provides a good level of purification 
and also results in a (usually) much smaller detector. 

Good energy resolution is required to prevent the tail of the (3(3{2v) spectrum from 
extending into the (3(3{Qv) region of interest. It also increases the signal-to-noise ratio, 
reducing the background in t he region of intere s t. Tw o-neutrino double-beta decay as 



background was analyzed by 



Elhott and Vogell (120021 ) 



• Ease of operation is required because these experiments usually operate in remote 
locations and for extended periods. 

• A large Qpp usually leads to a fast /3/?(0i/) rate and also places the region of interest 
above many potential backgrounds. 

• A relatively slow I3(3{2v) rate also helps control this background. 

• Identifying the daughter in coincidence with the (313 decay energy eliminates most 
potential backgrounds except /3/5(2z/). 

• Event reconstruction, providing kinematic data such as opening angles and individual 
electron energies, can reduce background. These data might also help distinguish light- 
and heavy-particle exchange if a statistical sample of l3(3{Qv) events is obtained. 

• Good spatial resolution and timing information can help reject background processes. 

• The nuclear theory is better understood in some isotopes than others. The interpre- 
tation of limits or signals might be easier for some isotopes. 

Historically, most j3j3 experiments have faced U and Th decay-chain isotopes as their 
limiting background component. A continuum spectrum arising from Compton-scattered 
7 rays, (3 rays (sometimes in coincidence with internal conversion electrons), and a parti- 
cles from the naturally-occurring decay chains, can overwhelm any hoped-for peak from the 
/3/5(Oz/) signal. This continuum is always present because U and Th are present as contami- 
nants in all materials. The level of contamination, however, varies from material to material 



35 



and even within different production batches of the same material. Over the years, great 
progress has been made identifying the level of contamination in materials that compose 
detectors and their shields. Once the location of the background-causing contaminant was 
determined, either cleaner materials were substituted or purification techniques were devel- 
oped. As a result detectors are now fabricated from amazingly pure components, some with 
activities as low as 1 /iBq/kg or less. 

As experiments reach for the {mpp) scale implied by the atmospheric neutrino 

oscillations, the rate of the signal is anticipated to be a few counts/ (t y). To observe such 
a small peak superimposed on a continuum of background, the background underlying the 
peak will need to be 1 count/ (t y) or less. This will require unprecedented low levels of U and 
Th within the experimental apparatus. In fact, the required levels of assay sensitivity are 
beyond what is currently technically achievable. Research and development are proceeding 
to improve assay capability in, for example, direct 7-ray counting, surface a- and /3-ray 
counting, and mass spectroscopy. 

At such extreme low levels of activity from U and Th, other components of the contin- 
uum that were previously a secondary concern may contribute a significant fraction to the 
background continuum. In fact it may be that no single component will dominate and a 
number of contributors may form a complicated mix. In the remainder of this subsection, 
we discuss the possible contributors and strategies to mitigate their effects. 

Double-beta-decay experiments are conducted deep underground to avoid cosmic ray 
interactions. At these depths, muons are the only surviving cosmic ray particles, but their 
interactions can produce high-energy secondaries of neutrons, bremsstralung 7 rays and 
electromagnetic showers. The experimental apparatus itself is usually contained within 
a detection system for muons. Hence any signal produced by the muon itself or any of 
it's prompt emissions from interactions in the apparatus will be eliminated by an anti- 
coincidence requirement. 

Neutrons, being neutral, can have sizable penetrating power through a shield. They 
generate background through (n, n'7) and (n, 7) reactions or (n, x) reactions that lead to 
radioactive nuclei. In particular background originating from (n, n'7) reactions within the 
experimental apparatus produce a rather non-descript spectrum. Any individual level in a 
nucleus might be excited only weakly, but there are many possible levels. Hence a significant 
7-ray flux without strong identiflable lines can appear, making, neutron induced backgrounds 



36 



difficult to diagnose. Neutrons have two sources: fission and (a,n) reactions in the cavity 
rock where U and Th levels can be relatively high, and muon interactions in the rock. Fission 
and {a,n) neutrons have energies less than about 10 MeV. These low-energy neutrons can 
be moderated and shielded from the detector with layers of hydrogenous material. The 
neutrons arising from muon interactions can have very high energy (>1 GeV) and therefore 
may penetrate the shield and induce a background-causing reaction near the detector. One 
particularly dangerous reaction is (n,n''y ) in P b; it can produce 7 rays with an energy 
very near to the of '^^Ge (IMei et all 120071 ). Even more onerous, is the Pb(n, n'7)- 



produced 3062- keV 7 ray, which has a double-escape peak very near Q/j/j. This single-site 
energy deposit could be problematic for Ge detector experiments using Pb as a shield. 
(See Section rVII.G.3l for a discussion of the Ge experiments.) Neutrons from muon-induced 
reactions can be reduced by going deeper underground. 

Fast neutrons, along with other hadrons, can produce radioac t ive nu clei that may create 



background for f3(3{0i') ( lAvignone et aU . I1992I : iBrodzinski et all 11990). For example, ^^Ge 
and ^°Co are produced in Ge and ^°Co is also produced in Cu. The production rates 
are significant on the surface of the Earth where the hadronic cosmic-ray flux is large. 
Because these nuclei have long half lives (0.8 and 5.3 y respectively), they will remain 
in the apparatus and potentially create background after the materials have been placed 
underground. These backgrounds can be mitigated by storing material underground and 
letting the isotopes decay away (not practical for Co) or by purifying the radioactive isotope 
from the host material and then minimizing its exposure above ground. Both strategies will 
be used in future efforts. Anthropogenic isotopes, suc h as ^°''Bi, have also been observed in 



jSP experiments (IKlapdor-Kleingrothaus et al. 



2004al ). Some exotic backgrounds have also 



been considered. For example, solar neutrinos are unlikely to be a significant background 
for most experimental configurations (lElliott and Engell . 12004 ) . 



B. Facility Requirements 

There is no doubt that a laboratory at significant depth is required for 13(3 experiments. 
Just how deep is a difficult question to answer, but clearly, deeper is better. However, lo- 
gistical concerns may induce a collaboration to choose a particular laboratory. The depth 
requirement for (3(3 is set by the flux of high-energy, muon-induced neutrons and subsequent 
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reactions that can lead to background. (See Section IVI.AI ) Depth -requirement calcula- 



tions are exceedingly difficult to perform. However, two recent works (IMei and Himd . 12006 ; 



Pandola et al. 



20071 ) have estimated the backgrounds in Ge-detector /?/?(0i/) experiments. 



Both found that the backgro und rate at a depth eq uivalent to that of Gran Sasso (3.1 km.w.e. 



flat overburden equivalent (IMei and Hime 



20061 )) could be kept below ~l/(t y), but that 



significant veto efforts would be required. The problematic fast neutron flux decreases with 
the muon flux, which, in t urn d rops by a factor of 10 for every 1.5 km.w.e. in depth. Mei 



and Hime (IMei and Hime 



20061 ) estimate that a depth of 5 km.w.e. (fiat overburden) would 



greatly reduce veto-system requirements. 

As the background requirements become more stringent, increasingly heroic efforts will 
be required to minimize activity internal to detector materials. In addition, as the experi- 
ments get larger and more complicated, the required infrastructure at the experimental site 
will become more demanding. For example, underground manufacturing of Cu to avoid 
cosmogenic activation of ^°Co, clean room facilities for assembly of sensitive parts, strong 
exclusion of Rn to prevent daughters from plating-out, and material purification capability 
are a few of the high-technology requirements for a remote underground location. Since 
many experiments (in addition to those on (3(3) will need an underground site and similar 
capabilities, the importance of a major facility to provide common infrastructure is clear. 
The opportunitie s and requirements for a deep underground laboratory are described in 

mm . 



Beier et al. 



C. Measurements to Constrain Nuclear Matrix Elements 



Part of the (3(3 program is an effort to reduce uncertainty in the nuclear-physics factors 
that, along with (mftFi), determine the rate of /3/5(Oz/). What data will be most useful? A 
recent workshop (jZuberl . l2005l ) addressed this question and summarized measurements that 
might help clear up nuclear-physics issues. Such measurments, says the report, should focus 
on high Qpp isotopes (the decay rate is proportional to Qf3/3^), even though low Q/?/? isotopes 
such as ^^'^Te and ^ss-q ^^.^ some interest. The 11 isotopes with a Qpjs greater than 2.0 
MeV are: ^^Ca, ^^Ge, ^^Se, ^^Zr, ^^^Mo, ^^^Pd, ^^^Cd, ^^^Sn, ^^^Te, ^^^Xe, and ^^"Nd. Most 
of these isotopes are currently part of a (3(3 experimental program. 

The Qj3[3 for these nuclei are known to only a few keV, with the exception of that in 
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^^Ge, which is known to ~ 50 meV. Precision measurements of the other Q- values would be 
useful because energy resolution is improving. The energy resolutions of the COBRA (~55 



keV) and CUO RE f^5 keV ) 



experiments are larger tha n or comparable to the uncertainty 



2003 



Redshaw et al. 



20071 ) for their chosen isotopes (4 keV for 



in the Q- values (lAudi et al. . 
^^^Cd and 2 keV for ^^"Te). It would also be helpful to measure the masses of certain nuclei 
involved in radiative EC-EC decay, because that decay rate can b e resonantly enhanced if 
the e nergy of the 2P-1S atomic transition is similar to the Q- value (ISujkowski and Wycechl . 



2004) ). At present the masses of the interesting nuclei are not known precisely enough to 



determine the degree of enhancement. 

Precise [3(3{2v) data in all nuclei used for /3/3(Oz/) and single-/3 data for the associated 
intermediate nuclei, where possible, are also needed. Rates for both process es have been 



used to fix the para r neter that plays a crucial role in QRPA equations (IRodin et al. 



2003 



2006 



Suhonen. 



20051 ). (See Sec. IV.BI ) More data would test internal consistency, and 



help decide which Qpp calibrator is better. 

In (3[3{2v), and in [3[3{^v) as well in the closure approximation is not used, the expressions 
for the overall matrix element (see Eqs.f l38l) and ( 124|) ) contain matrix elements of the charge 
changing weak current from both the parent and daughter nuclei to states in the interme- 
diate nucleus. In the long-wavelength approximation, which the high-momentum virtual 
neutrino prevents for (3(3{Qv), these reduce to the Fermi and Gamow- Teller operators. GT 



strength distributions from the parent nuc 



with charge-exchange reactions (lAmos et al. 



eus in the releva nt isotopes can be measured 



2007 



Ejiri 



2000l ) , and should be measured for 



all the 9 high-Q- value isotop es of experimental interest. Extracting GT strengths from the 
daughter (lAmos et all 120071 ) and the strengths for higher-multipole analogs of the GT oper- 
ator is critical for /?/3(0z/) but much harder because the relations between the cross sections 
and operator matrix elements are more complicated, and polarized beams are sometimes 
required. Transitions from the parent can be studied, e.g., with (p,n) or (^He,t) reactions 
and transitions from the daughter with (n,p), (d,^He), or (t,^He). But even if the connec- 
tion between cross sections and transition strengths can be properly made, the strength is 
still the square of a matrix element, so that the corresponding sign, important for (3(3, is 
unavailable. 

Muon capture is governed by matrix elements similar to those connecting the intermediate 
and daughter nuclei in (3(3. Because the muon is heavy, states with any multipole can be 
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populated during capture, and partial muon captur e rates to specific states can be l in ked to 



2002h 



Volpe 



the co rresponding parts of the matrix elements (iKortelainen and Suhonen 
( I2OO5I ) proposed charged- current neutrino-nucleus scattering as a technique to constrain 
Mqi/. By using both v and 17, one can probe matrix elements of one-body operators to 
intermediate states from both the parent and daughter nucleus. 

The study of pair correlations of nucleons in the ground state of the parent and daughter 
nuclei allows one to probe aspects of the wavefunctions that are important for (5(5 matrix 
elements even when the closure approximation is used. Double-beta decay can be represented 
as the removal of two neutrons followed by the addition of two protons, and calculations 
show that when the representation is decomposed into multipoles, the J'^ = 0"^ channel is 



the m ost important. Recent results of (p,t) reactions for the Ge-Se system (IFreeman et al. 



20071 ) have shown no evidence of pairing vibrations that have been found in low-A systems. 
^^Ge and ^^Se have quantitatively similar neutron pairing correlations, and calculations will 
be easier there than they would be otherwise. 

Other processes, such as pion double-charge exchange, electromagnetic transitions 



to i sobaric analogue stat e s (lEiiril. 



ter (lAunola and Suhone 



2 



1996 : 



2006ar). and (5l5{2u) t o exc ited state in the daugh- 



Suhonen and Civitaresd . 



19981 ) can also help constrain 



/3/3(Oi/) matrix elements. But despite the usefulness of all these auxiliary measurements, 
none will provide a magic bullet. Such data will help judge the quality of calculations, but 
not reduce uncertainty to near zero. For that, the best data are probably Tj^g themselves, 
the same quantities one is trying to calculate. If the lifetimes are measured in several nu- 
clei, one can check directly whether the calculations work, most simply by fixing 
from one nucleus and seeing whether the other lifetimes are correctly predicted. Of course, 
this assumes that light-neutrino exchange rather than some other mechanism is driving the 
decay. 



VII. FUTURE PROGRAM 



A. Previous Experiments 



Since the first direct observation of (5(5{2v) in 1987 ( lEUiott et ali 119871 ) . decay rates 



in a number of other isotopes have been measured, either directly or with geochemi- 
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cal/radiochemical techniques. The recent work of Rodin et al. ( 



Rodin et al. 



2OO3I . l2006h 



has rekindled interest in the precise measurements of T'lj2- iBarabashI (j2006a|), in a review of 
the /3/3(2z/) experiments, has recommended average values, which we quote in Table [ITl 
Here we just note that because is proportional to Qjsp^^, whereas Vqi, is proportional 
to (5/3/3^, T^y might be larger than r2,yif Qpp is small and geochemical measurements of 
the decay rates that do not distinguish between /3/3(Oz/) and (3[3{2i') might still produce 
a competitive limit on (rrifif ^). This is the case in ^^ ^Te, where the given in Table HTl 
produces a limit on {mpp) (IBernatowicz et ali Il993l ) that is onl y slightly worse than the 



best limits given in Table IIIII The nucleus ^^^U is another case (ITurkevich et al. 



199l|) in 



which the radiochemical experiment could not distinguish f3(3{Qv) and (3(3{2v). ^ss^j (iecays 
to ^'^'^Pu, which in turn a-decays with a 87.7-y half life. A sample of U salt had been stored 
for 33 years and was milked for its Pu content. By counting the ^^'^Pu a, the (3(3 half 
life was determined. Because of its low Qpp the high observed decay rate (compared to the 
theoretically predicted TIJ2) has been interpreted as evidence for (3(3{Qi'). Unfortunately, the 
experiment is difficult to repeat because of the special nature of the U-salt sample available 
to the experimenters. 

The entire history of (3(3 mea surements up to about 2001, includin g /3/5(2z/), /5/5(Oz/), and 
Majoron modes, can be found in lTretyak and Yu G. Zdesenkd (120021). The best experim ents 



to date have been in Ge. The Heide lberg-Moscow (iKlapdor-Kleingrothaus et all l2001al ) and 



IGEX f lAalseth et al 



2002a 



20041 ) experiments have provided the best limits on Tf^ by 



building detectors with the lowest backgrounds. 



B. Overview of the Future Program 

In the inverted hierarchy with a mughtest near meV, {mjsp) will be near 20-50 meV. 
The Tfj2 resulting from this neutrino mass will be near 10^'' y with a resulting count rate of 
few/(t y). On the other hand, if the recent claim of {mpp) ~ 400 meV is borne out, the half 
life will be nearer to 10^^ y, with a count rate of a few 100/ (t y). A precision measurement 
(20% or better) requires a few hundred events. To make such a measurement in a 2-3 year 
run, it seems prudent to build an experiment containing 100-200 kg of isotope that can 
be expanded to 1 ton at a later time. But collecting statistics in a reduced background 
experiment is not the only goal. The recent claim for /3/3(Oz/) is controversial. Future claims 
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must provide supporting evidence to strengthen the argument (lEUiottl . l2006l ) . This evidence 
might include any number of the following. 



• To show that /5/3(Oz/) likely exists, one needs a combination of: 

— a clear peak at the correct PP^Ou) energy, 

— a demonstration that the event is a single-site energy deposit, 

— measured event distributions (spatial, temporal) representative of /3/3(Oi^), 

— a demonstration that the measured decay rate scales with isotope fraction. 

• To present a convincing case, one needs: 

— an observation of the 2-electron nature of the /3/3(0z/) event, 

— a demonstration that the kinematic distributions (electron energy sharing, open- 
ing angle) match those of /3/3(0z/), 

— to observe the daughter nucleus together in time with the /3/3(Oi/) decay, 

— to observe excited-state decay with parameters indicating /3/3(Oi/). 

• To remove all doubt, many of the above (3(3{0i>) indicators should be measured in 
several isotopes. 

The proposals summarized below in Section IVII.CI address these possibilities in different 
ways. 

1. The Number of Required Experiments and Their Precision 

If /?/3(0i/) is observed, the qualitative physics conclusion that neutrinos have a Majorana 
nature will have a profound effect on the development of models for mass generation. Quan- 
tifying the physics associated with lepton-number violation will be more challenging because 
of the uncertainty in Mq^. If the matrix elements were precisely known, one could compare 
Foi/ measurements in several nuclei to address the question of the underlying mechanism. 
If, in addition, the mechanism was certain, one could determine the LNVP from a lone Fot/ 
measurement. 
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Recently a number o f papers have appeared ( 



Deppisch and Pad . 



2007 



Gehman and Elliott 



Bahcall et al. 



2004a 



Bilenky et al. 



2004; 



20071 ) that try to quantify, in light of matrix- 



element uncertainties, the number and type of experimental results needed to draw quanti- 
tative conclusions. Although the papers make different assumptions, the concensus is that 
if (3j3{Qv) exists, measurements of Fg^y in at least 3 different isotopes are warranted. All 
the analyses assume that existing calculations provide an estimate of the difference in rates 
produced by different mechanisms, and require that the experiments determine the rate 
to within 20%. This last number comes from the level of variation in calculated matrix 
elements. 



2. Kinematic Distributions 



If /3/3(Oz/) is observed, future experimenters will want to study the k inematic distr i butions 
associated with the outgoing electrons. The NEMO 3 experiment ( Arnold et al. . boOSa ) 
has shown the power of this type of analysis with (3(3{2v). The collaboration used the 
spectrum of electron spectra from f3(3{2v)'m ^'^°Mo, taken one at a time (the lone electron 
spectrum)., to show that the decay proceeds predominately through the 1^ ground state of 
the inter mediate ^°°Tc nucleu s. The effect on the spetrum of this "single-state dominance" 
is small (jSimkovic et a/.l . 120011 ) . so 10^ events were required for the analysis. But the analysis 
demonstrates that important information contained within kinematic distributions. 

Figure [6] shows an example of the kinematic distributions for /3/3(0i/). In this Figure, 
the lone electron spectra are shown for two possible exchange mechanisms. The mountain 
shaped curve is for light-neutrino exchange mechanism and the valle y shaped curve is for 
right-handed interactions in both the leptonic and hadronic currents ( iDoi et all Il985l ). It 
is clear from this Figure, that only a few 10s of events would be needed to show that one of 
these mechanisms dominates the decay. However, a large sample of events would permit an 
analysis to constrain the fraction of a competing mechanism. The opening-angle distribution 
is another observable that is sensitive to aspects of the underlying exchange mechanism. 

The only detectors able to regist er these types of distribu tions involve thin-foi l, tracking- 



volum e sandwiches. The MOON (iNakamura et alx . l2006al ) and SuperNEMO (IBarabashl . 



20041 ) proposals are examples. There are competing demands in the design of such exper- 
iments, one would like a thick foil for a large source mass, but a thin foil to minimize 
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energy loss and scattering so the integrity of the distributions will be retained. As a result, 
such experiments will have to be large with many electronic channels in order to collect the 
required numbers of events. 



C. The Experiments 

In this section, we describe a number of proposed experimental programs for /5/?(0z/). A 
summary of the proposals is given in Table llVi 

1. CUORE (Cryogenic Underground Observatory for Rare Events) 



Before addressing the CUORE experiment, it is appropriate to discuss the general prin- 
ciples of cryogenic, or thermal, detectors, and the pilot experiment, CUORICINO. That 
detector is an array of bolometers made from a material containing the parent decay iso- 
tope. In CUORICINO and CUORE the individual bolometers are 5x5x5 cm^ single 
crystals of Te02. They are made from natural- abundance Te, which is 33.8% ^^°Te. They 
are operated at very low temperatures where they have tiny specific heats. They are weakly 
thermally coupled to copper support frames that act as the heat bath. When an event 
leaves energy in the crystal, the temperature rises and that rise is measured by a sensitive 
thermistor, which produces a signal proportional to the energy deposited. The heat leaks 
back through the weak thermal couplings and the signal decays as the temperature returns 
to its baseline. 



This technique was suggested for j3(3 decay searches bv lFiorini and Niinikoskil (1984r), and 



applied earlier by th e Milano group in the M IBETA experiment (lArnaboldi et al 



fore CUORICINO flArnaboldi et al 



200 



20031) be- 



The CUORICINO bolometers are dielectric 



and diamagnetic, and are operated at temperatures between 8 and 10 mK (lArnaboldi et al. 



2003 



2004J . l2005l ) . According to the Debye Law, the specific heat of Te02 crystals is given 
by C{T) = P{T/Qd)^, where /3=1994 J/(K mol) and Qd is the Debye temperature. In 
these materials, C(T) is due almost exclusively to lattice degrees of freedom. The con- 
stant 0£) was specially me a sured by the Milan group for 5x5x5 cm^ Te02 crystals 



as 232 K fArnabo 



272K flWhite et al. 



di et al. 



20031 ). which differs from the previously published value of 



19901 ). The specific heat follows the Debye Law down to 60 mK. The heat 
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capacity of the crystals, extrapolated to 10 mK, is 2.3x10"^ J/K. With these values of the pa- 
rameters, an energy deposition of a few keV will result in a measurable temperature increase 
AT. The temperature increase caused by the deposition of energy equal to the total /3/3(0i/) 
decay energy, Qp/s = 2530 keV, would be 0.177 mK. To obtain usable signals for such small 
temperature changes, very sensitive thermistors are required. In CUORICINO, AT is mea- 
sured by high resistance germanium thermisters glued to each crystal. More details can be 



found in 



1997 



Fiorini and Niinikoskil (jl984j ). and in subsequent publications (lAUessandrello et al. 



19981). 



The thermistors are heavily doped high-resistance germanium semiconductors, with an 
impurity concentration slightly below the metal-insulator transition. High-quality thermis- 
tors require a very homogeneous doping concentration. For the CUORICINO thermistors 
(Neutron Transmutation Doped (NTD)), this was achieved by uniform thermal neutron 
irradiation throughout the entire volume in a nuclear reactor. The electrical conductiv- 
ity of the devices depends very sensitively on the temperature because they use variable 
range hopping (VRH) mechanisms. The resistivity varies with temperature according to 
p = poGxp (Tq/T)'^, where the constants po, Td and 7 all depend on the doping concentra- 
tion. In the case of thermistors operating with VRH mechanisms, 7 = 1/2. 

Thermistors can be characterized by their sensitivity, A{T), defined as follows: A{T) = 
\d{lnR) / d{lnT)\ = •yiTo/Ty, where the resistance is -R(T) = Rq exp (To/Ty. The parameter 
Rq = Po{d/a), where d and a are the distance between the contacts and the cross section 
of the thermistor, respectively. The values of Ro, Tq and 7 must be measured for each 
thermistor. This is done by coupling the thermistor to a low-temperature heat sink with a 
high heat-conductivity epoxy. The base temperature of the heat sink is between 15 and 50 
mK. A current flows through the device and a V-I load curve is plotted. The curve becomes 
non linear because of the power dissipation and dynamic resistance, causing the slope of 
the I(V) curve to change from positive to negative. The optimum operating biasing current 
occurs where dl/dV increases rapidly with increasing bias voltage, Vb, maximizing the signal 
to noise ratio. The parameters of each thermister are determined from a combined fit to a 
set of load curves at different base temperatures. The charact erization process was described 



in deatil earlier for Si thermistors (AUessandrello et al. 



19991 ) ■ and the same process is used 



for the CUORICINO Ge thermistors. 

The thermistors in CUORE will be the same as those used in CUORICINO; they were 
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produced by the UC Berkeley/LBNL group (IHaller et g/J . Il982l ). It is necessary to optimize 
the neutron doping of the Ge. This is done by using foils of metal with long-lived (n, 7) 
radioactive daughter nuclides irradiated along with the Ge to calibrate the neutron flux. 
Accordingly, the neutron exposure can be accurately evaluated without having to wait for 
the intense radiation of the ^^Ge in the sample to decay. Following the decay period, the Ge 
is heat treated to repair the crystal structure, and cut into 3x3x1 mm strips. Electrical 
connections are made by two 50 fim gold wires, ball bonded to metalized surfaces on the 
thermistor. The thermistors are glued to each bolometer with 9 spots of epoxy, deposited 
by an array of pins. 

Figure [7] shows the CUORICINO structure. Each of the upper 10 planes and the lowest 
one consists of four 5 x 5 x 5-cm'^ crystals of natural isotopic abundance, as shown in the 
upper right hand figure, while the 11th and 12th planes have 9 3 x 3 x 6-cm^ crystals, as 
shown in the lower right hand figure. In the 3x3x6 cm^ planes the central crystal is fully 
surrounded by the others. 

The crystals for CUORE will be prepared the same way as for CUORICINO. The crystals 
were grown with pre-tested low-radioactivity material by the Shanghai Institute of Ceramics 
of the Academy of Science (SICAS) and shipped to Italy by sea to minimize activation 
by cosmic rays. They were lapped with a specially selected low contamination polishing 
compound. All these operations, as well as the mounting of the tower, were carried out 
in a nitrogen atmosphere glove box in a clean room. The mechanical structure is made of 
OFHC Copper and Teflon, both of which were previously tested to be sure that measurable 
radioactive contaminations were minimal. Thermal pulses are recorded by the NTD Ge 
thermistors thermally coupled to each crystal. The gain of each bolometer is calibrated and 
stabilized by a 50-100 kQ resistor attached to each abs orber and acting as a heater. Heat 
pulses are periodically supplied by a calibrated pulser (lArnaboldi et a/.l . |2004| ). The tower 
is suspended from the 50 mK plate of the cryostat through a 25 mm copper bar attached to 
a steel spring. The steel spring provides mechanical isolation from vibration of the cryostat 
that can result in heating and spurious pulses in the detector. 

The CUORE detector will be an array of 19 towers, each similar to the CUORICINO 
tower except that all crystals will be 5 x 5 x 5 cm^ (See Fig. [8]). It will have 988, 760-g 
bolometers, containing a total of about 750 kg of Te02, or ^200 kg of ^'^°Te. The expected 
background is projected from the CUORICINO background, and the improvements thus 
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far achieved by the collaboration, as well as the effect of coincidence cancellation due to 
the granularity of the detector. The predicted background rate for the first phase of the 
experiment is 0.01 counts/(keV kg y). Substituting into Eq. with n„ = 1, a = 0.348, 
e = 0.84, W = 162, M = 760 kg, t = 10 y, 6{E) = 7 keV and b = 0.01 counts/(keV 
kg y) the predicted sensitivity is T{y2(^'^°Te)^ 2.5 x 10 ^^^ y. According t o three recent 
nuclear structure ca 
10"^^/v JPoves et al. 



culations, Fn = (4.84^^;^^) x lO^^V y (IRodin et al. 
20061 ). or = 5.13 x 10"^^/y (ICivitarese and Suhoned . 



2006), Fm = 2.57 X 



2003 



These results and the sensitivity above yield (m^^)-limits of 47 meV, 53 meV, and 45 meV, 
respectively. A version of CUORE that is isotopically enriched to 80% in ^^'^Te would reduce 
these values by a factor of 1.54 to 31 meV, 34 meV, and 29 meV, respectively, which covers 
the inverted hierarchy if mughtest is meV. The aim of the CUORE collaboration is to reduce 
the background further. 

The CUORE experiment is approved and funded by the Instituto Nazionale de Fisica 
Nucleare (INFN) in Italy, and was recently approved for construction and put in the 2008 
fiscal-year budget of the US Department of Energy. However construction has already begun 
in Laboratori Nazionali del Gran Sasso (LNGS), in Assergi, Italy with Italian funds. Data 



Ardito et al. 



acquis ition is scheduled to begin in 2011. A detailed description is given in 
(bood). The recent resuhs from CUORICINO with a total exposure of 11.8 kg-y of ^^^Te 



provide a lower limit of TP^ ^ 3.0 x 10^^ y (lArnaboldi et al. 



20071 ). Using the nuclear 



-1/2 

structure factors listed above, the bounds on {mpp) are 0.42, 0.58, and 0.40 eV, respectively. 
The mission of the first phase of CUORE, then, is to improve the half-life sensitivity by a 
factor of 100, and that of {mpj^) by one order of magnitude. 



2. EXO (The Enriched Xenon Observatory) 



EXO is a proposed large next generation experiment that will search for the f3(3{Qv) in 
^^^Xe. The goal is to use between one and ten tons of Xe enriched to 80% in ^^^Xe in a 
tracking Time-Projection Chamber (TPC). The intent of the collaboration is to trap and 
identify the daughter ^^^Ba ion by laser spectroscopy. The decay results in a ^^^Ba+"^ ion 
that will rapidly capture an electron, leading to leading to ^^^Ba+ ion which is stable in Xe. 
These ions can be positively identified via their atomic s pectroscopy by optical pumping 
with blue and red lasers. This technique was suggested by iMod (1l99ll ) and is described in 
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the context of this experiment by Danilov et al. (iDanilov et ali l2000al lbl). To exploit this 
method of positively identifying the daughter nucleus, it will be necessary to capture the 
ion at the decay site and transfer it to an ion trap by optical pumping. Alternatively, the 
ion could possibly be identified at the decay site by directing the lasers. If either of these 
techniques works, backgrounds that do not result in the correct ion will be eliminated. The 
development of these techniques is a significant research and development challenge, but one 
with a very high potential payoff. 

The atomic level structure of the Ba ion is shown in Fig. [91 The 493 nm ground-state 
transition is very strong allowing the ions in the 6^S'i/2 ground-state to be optically excited 
to the 6^Pi/2 excited state. The decay of this state has a 30% branching ratio to the 5^/^3/2 
metastable state. The Ba'*' ion is identified by irradiating this state with 650 nm red 

laser light, exciting the ion back to the state, which decays with a 70% branching 

ratio to the ground state by emitting 493 nm blue light. The excited state has a mean 

life of 8 ns and when saturated can emit ~10^ 493 -nm photons per second. This wavelength 
is compatible with the maximum quantum efficiency of bialkali photomultiplier tubes. The 
bright glowing spot shown in Fig. [9] is a photograph of an actual excitation in a test trap. 

The isotopic enrichment of tons of xenon is technically feasible because xenon is a noble 
gas that can be directly introduced into a mass-separating centrifuge with no chemistry. 
In addition, the residual xenon can be returned for normal industrial use. This represents 
a significant cost saving in the acquisition of the raw xenon. The xenon gas can also be 
continuously purified during the operation of the experiment, in particular to eliminate the 
ubiquitous radioactive radon and krypton isotopes. 

There are two activities that will presumably lead to the final construction of the ton-level 
EXO with barium-ion tagging. One is the development of the tagging procedure discussed 
above, and the other is the development of a high-resolution liquid xenon TPC with good 
tracking ability. The EXO collaboration has approval and funding to construct a 200-kg 
hquid TPC with Xe enriched to 80% in ^^exe. This project is called EXO-200. It will not 
have the barium-ion tagging feature; however, it is intended as an R&D program for the TPC 
itself. The enriched isotope has been purchased and is in house. This construction of EXO- 
200 is presently underway, and the experiment will be installed and operated in the Waste 
Isolation Pilot Project (WIPP) site in Carlsbad, New Mexico. This site has an overburden 
of 1600 m.w.e. to reduce the backgrounds associated with cosmic rays. The detector will 
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have an estimated 115 kg of fiducial volume contained in a thin- walled copper vessel 40 
cm in length. The copper was supplied by a German company, and has been assayed by 
Inductively Coupled Plasma Mass Spectroscopy (ICPMS). It contains the following limiting 
levels of radioactivity: "^°K < 6 ppm, ^'^^Th < 0.5 ppt and ^ss^j < q 3 ^jj-^g copper for 
the TPC was stored at the Deutches Electron Synchrotron Laboratory (DESY) in Hamburg, 
Germany, in a concrete bunker to reduce activation from energetic neutrons produced by 
cosmic ray muons. It has now been shipped to Stanford. It is possible that in the future, 
the collaboration will use Teflon for the vessel; however, the first vessel will be made from 
this copper. An artist's conception of the EXO-200 is shown in Fig. [101 

The high voltage cathode will be located in the middle of the TPC volume. The design 
calls for up to 3.5 V/cm. At each end there are 114 x-wires and 114 y-wires (actually at 
60° pitch with respect to the x-wires). These are for charge collection, energy, and position 
information necessary for the reconstruction of the trajectories of the electrons. A 1-MeV 
j3 particle creates about 50000 electrons in the liquid xenon. One interesting property of 
xenon gas and liquid is that they are efficient scintillators. The scintillation light will be 
detected by 258 large area photodiodes on each side of the TPC. The use of both collected 
charge and scintillator light improves the energy resolution. The combination was tested in 
a small chamber, and the data imply that EXO-200 will have an energy resolution of about 
1.5%. 

The installation in the WIPP site is scheduled to occur in mid 2007. The cryostat has 
been constructed and tested at Stanford University at the time of this writing. The EXO 
Collaboration anticipates reaching a sensitivity of about 0.3 eV for {mpp). The actual design 
of the full-size EXO detector will occur after the research and development of the barium- 
tagging scheme are completed. 



3. Next Generation "^^ Ge Double-Beta Decay Experiments: Majorana and GERDA 



As st ated earlier, the most sensitive (5[3 



Moscow (IKlapdor-Kleingrothaus et al. 



Oi/) experiments so far are the Heid e 



2001a 



and the IGEX (lAalseth et al 



2002 



berg 



20041 ) 



experiments. In addition, the claimed observation discussed earlier was ba s ed on a reana l- 



ysis of the Heidelberg-Moscow data (IKlapdor-Kleingrothaus et al 



Klapdor-Kleingrothaus and Krivosheinal . 



2001al . 



2003 



2004a 



20061). The goals of the Majorana and GERDA 
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^^Ge experiments are twofold, first to confirm or refute the claimed evidence, and second, 
in the case of a positive result, to make a measurement of the half life with a 20% or better 
uncertainty. If the claim is not correct, the goal is to make the most sensitive search for 
this decay. We will discuss the Majorana experiment first because it will use improved 
ultra-low background cryostat technology similar to that used in the IGEX experiment, 
which, except for the copper, was similar to that used in the Heidelberg-Moscow experi- 
ment. By contrast, the GERDA Collaboration is pursuing a novel technique of cooling and 
shielding the bare Ge detectors by directly immersing them in liquid argon (IHeusserl . Il995l ) 
and water as low-Z shielding. This technique was the cornerstone of a proposal by the GE- 



NIUS (GErmanium Nitrogen Undergroun d Setup) collaboration (iKlapdor-Kleingrothaud . 



1997 



Klapdor-Kleingrothaus et al 



Sasso in the GENIUS test facility flKlapdor-Kleingrothaus et al 



200 Ibl) and was tested in ari initial experiment at Gran 



2003). Simi 



Zdesenko et al. 



arly the GEM 



200ll ) also involved 



(Germanium Experiment for neutrino Mass) proposal (ji 
bare Ge detectors, but with pure water forming the shield instead of liquid nitrogen. Neither 
the GEM nor GENIUS collaboration is currently active. In the end, the Majorana and 
GERDA collaborations envision a joint one-ton experiment, using whichever of their two 
experimental techniques proves to be the most appropriate to reach the final goals. 

^^Ge decay has several significant advantages. First, all three recent nuclear structure 
calculations are in fair agree ment and predict large nu clear structure factors, - Fjy. They are : 

2Q06|k^ 



200.1 



Fn = 1.22l^:i^ X 10"^^ /y fiRodin et al 
20061 ) : and = 7.01 x 10^^^ y (jCivitarese and Suhonen 



4.29 X 10-^Vy flPoves et al 



2003al jbl). Fat is proportional 



to the square of the nuclear matrix element, so the matrix element differences themselves 
are not large, whereas in previous estimates they varied by a factor of 3. It is generally 
believed that more reliable shell model calculations will be possible in the future for this 
nucleus. In addition, ^^Ge experiments have the lowest background demonstrated thus far. 
The background rates at IGEX and Heidelberg- Moscow were ^0.06 counts/(keV kg y), after 
the application of pulse-shape analysis. Many of the backgrounds have been identified, and 
both future experiments are projecting significant reductions. Finally, Ge detectors have 
the best energy resolution of any of the proposed experimental techniques for double-beta 
decay experiments. Resolutions of 3 keV at the decay energy of 2039 keV, or ~0.15% full 
width at half maximum (FWHM), have been achieved. The main difference between the 
Majorana and GERDA proposals are the method of cooling the detectors, and shielding 
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and background control. Both propose to use detector segmentation and pulse shape dis- 
crimination to reduce backgrounds. The goals of both proposed experiments is to probe the 
quasi-degenerate neutrino mass region above 100 milli-electron (meV) region in the initial 
phases. In a subsequent phase, and possibly with collaboration between Majorana and 
Gerda, the goal is a one-ton experiment to explore the inverted hierarchy down to about 20 
meV. 

To achieve these sensitivities it is necessary to fabricate the detectors from germanium en- 
riched to 86% in ^^Ge. This material can be produced in Zelenogorsk, Russia, in the required 
quantities via centrifuge. The Ge metal is converted into a GeF4 gas that is very stable at 
high temperatures, but must be converted into an oxide or metal after centrifugation. This 
is the same type of starting material that was used in both the Heidelberg-Moscow and 
IGEX experiments. No observable internal contamination of the Ge was found in either of 
those experiments. 

The Proposed Majorana Experiment: The original Majorana concept was for 
eight modules with fifty-seven, 1 .1 kg segmented Ge d etectors each, for a total of ~500 kg 
of Ge, enriched to 86% in ^^Ge f Gaitskell et all . 12003 ). The detectors within each module 
would be arranged in a hexagonal configuration of 19 towers, each three detectors high, 
inside of a copper cryostat. The cryostats will be electroformed as in the case of the IGEX 
cryostats. An intense R&D project is underway to reach the following levels of radio-purity 
for the Cu: ^^^Bi: 0.3 /iBq/kg, and ^^^Th: 0.1 /iBq/kg. These two isotopes produce the 
most serious of backgrounds, for example the 2615-keV 7 ray from ^os^l at the end of the 
^^^Th chain. The background associated with this 7 ray is difficult to eliminate sufficiently 
by analysis as discussed below. For this reason, the reduction of thorium, as well as the 
^^^Bi, in the copper is a critical step. 

The electroforming of the copper cryostat parts begins with the selection of pure materi- 
als. The copper is electroplated on stainless steel mandrels from a copper sulfate solution. 
The solution is prepared with semiconductor-grade acids. The copper sulfate is purified by 
re-crystallization. This procedure was tested by artificially contaminating the copper sul- 
fate, and then re-crystallizing it several times, measuring the radioactivity at each step. The 
solution is constantly re-circulated through micro filters to remove oxides and precipitates. 
It is also circulated through a barium scavenge filter to remove radium. The plating tanks 
have a cover gas to reduce oxides. The copper parts are machined periodically to prevent the 
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formation of porous copper, which can lead to vacuum leaks and poor structural strength. 
While this is a time consuming process, the payoff is copper cryostat parts with less than 1 
/iBq/kg of 232Th. 

Each module of 57 detectors comprises 19 three-detector towers as shown in Fig. [TTl 
The final design of the individual Ge detectors has not yet been completed; however, the 
base plan has 1.1-kg segmented detectors with a geometry similar to that of the Canberra 
CLOVER detectors. Some initial studies concerning the jo int use of seg r aenta tion and pulse 



shape analysis with these detectors has been encouraging (lEUiott et al. 



20061 ). Background 



events from 56,57,60qq^ ^^Zn, and ^^Ge, formed in the Ge crystals by spallation reactions from 
high-energy neutrons produced by cosmic ray muons can be identified event-by-event by the 
deposition of energy in more than one detector segment, and by pulse-shape discrimination. 

Another significant source of background is externally produced 7 rays. Gamma rays 
that can deposit 2 to 3 MeV in a Ge detector most likely scatter several times before either 
being absorbed or leaving the crystal. This creates multiple sites where clouds of electrons 
and holes are formed. The result is a complex displacement current pulse that can be dis- 
tinguished from one that results from a single-site event like that expected from /?/5(0z/). 
Sample experimental pulses are shown in Fig. [T2j The figure displays an experimental spec- 
trum of background after segment cuts and pulse-shape discrimination have been applied. 
All of the 7 ray lines have been significantly reduced except for the 1593-keV double-escape- 
peak line of the 2615-keV 7 ray in ^os'pi This well-known line is from the single site creation 
of electron-positron pairs, in which both annihilation 7 rays completely escape the detec- 
tor. One also sees that about 1/3 of the continuum remains because it results from single 
Compton scatters. 

Accordingly, by lowering the background in the copper cryostat parts, by carefully screen- 
ing and selecting the front-end electronic parts for low background, and by applying seg- 
mentation and pulse-shape discrimination, the Majorana Collaboration hopes to reduce 
the background in the region of interest at 2039 keV by a factor of 150 below that of IGEX 
or Heidelberg-Moscow. The target half-life sensitivity, after an exposure of 0.46 ton- years. 



is 5.5 X 10^^ y. I 



y (IRodin et al. 



we use the central value of the nuclear structure factor, F/v = 1.22 x 10 / 
:etime corresponds to a (ttibfi) ^61 rn eV to a 90% C.L. If 



2003 



20061), the h 



the claimed observation is correct (iKlapdor-Kleingrothaus et a/.U2004bl ). then this exposure 



would result in a measurement of the half life with an accuracy of 10%. 
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Majorana is a collaboration of about 70 physicists from 16 institutions in Canada, 
Japan, Russia and the USA. At the time of this writing, it has been approved for major 
research and development. 

The GERDA Experiment (GERmanium Detector Array): The Germanium De- 
tector Array (GERDA) is a planned array of Ge detectors fabricated from germanium en- 
riched to 86% i n ''^Ge, cooled and shielde d by direct immersion in liquid argon or nitro- 



gen (lAbt et al 



2004 



Schonert et al. 



20051 ) . The facility is presently under construction in 



the Laboratori Nazionali del Gran Sasso (LNGS) in Assergi, Italy. It is planned in three 
phases. The first phase will proceed with the 5 detectors from the Heidelberg-Moscow (HM) 
Experiment, with 11.9 kg of enriched Ge, and the 3 detectors from the IGEX Experiment 
with 6 kg of enriched Ge. All of these detectors are closed-end p-type semi-coaxial detectors, 
and are unsegmented. In their copper cryostats, these detectors achieved background rates 
of :^0.06 counts/(keV kg y) after pulse shape discrimination that reduced the number of 7- 
ray events producing multi-site pulses. The detectors have been underground since the early 
1990's, the HM detectors in LNGS, and the IGEX detectors in the Laboratorio Subteranio 
de Canfranc in Spain. The goal of GERDA is to reduce the background to levels between 
0.01 and 0.001 counts/(keV kg y) in the initial phases via ultra-clean materials, pulse-shape 
analysis, segmentation, and live and bulk shielding with liquid argon. An artist's conception 
of GERDA is shown in Fig. [H 

Phase I, comprising an active mass of 17.9 kg, is supposed to achieve a half-life sensitivity 
of ^ 3 X 10^5 y at (90% CL). For Phase II, an additional 37.5 kg of Ge, enriched to 87% in 
^^Ge, has been delivered to the Max-Planck-Institute in Munich in the form of Ge02. It will 
be stored underground until fabrication begins. The collaboration is testing one 6-fold seg- 
me nted p-type dete ctor, and one 18-fold segmented n-type detector for possible use in Phase 
II ( lAbt et all 120071 ). The goal of Phase II is to reach a half-life sensitivity of 1.4 x 10^^ y 



(90% CL) and an exposure of 100 kg y. 
factor Fn = 1.22 x lO^^Vy Jpodin et al. 



f we use the central value of the nuclear structure 



2003 



20061). this corresponds to a (rriM) of ^ 124 



Klapdor-Kleingrothaus et al\ ( l2004bl ). is 



meV. In the case that the claimed observation by 
correct. Phase I would observe it and Phase II would make an accurate measurement of the 
lifetime. Fig. shows an artist's conception of a Phase II detector tower that would hang 
directly in the liquid argon. 

The GERDA cryostat is a vacuum-insulated stainless cryostat with an inner copper liner. 
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The outer cylinder is 4.2 meters in diameter and 8.9 meters long. The inner volume is 70 
m^, the mass of the empty vessel is 25 metric tons and it will contain 98 tons of liquid 
argon. There will be a clean room over the vessel and a rail system to lower and position 
the individual detector strings. 

Monte Carlo calculations were performed in collaboration with members of the Majo- 
RANA computation group to predict the background levels. GERDA will operate in Hall-A 
of the LNGS with an overburden of ~ 3600 m.w.e. Accordingly, it will have a muon veto 
shield to reduce backgrounds associated with cosmic-ray muons. Given the low-Z material 
used as shielding in the vicinity of t he germanium array, the muon induced background will 
be 0.0001 counts/(keV kg y) or less (jPandola et aU 120071 ). In addition to segmentation cuts, 
and pulse-shape analysis, anti-coincidence signals between individual Ge detectors can be 
used to tag and eliminate external 7 rays. Coincidences from 7 rays in decay chains can be 
used to further reduce backgrounds. Finally, R&D toward the use the scintillation light of 
liquid argon as an activ e veto system is underw ay. The technique may be implemented in a 
later phase of GERDA flDi Marco et all kom . 

The ultimate goal of a world wide collaboration would be a half-life sensitivity of 10^*^ y, 
corresponding to ~25 meV. This would presumably be an experiment with the order of one 
ton of enriched germanium. This endeavor may very well be a collaboration between the 
GERDA and Majorana groups, depending on which technology proves to be the best with 
respect to background, detector stability, etc.. There exists a memorandum of understanding 
to that effect, and it is clear that there will very probably be only one Ge experiment in the 
world of that magnitude. 

GERDA is a collaboration of 80 physicists from 13 institutions in Belgium, Germany, 
Italy, Poland, and Russia. It was approved in November, 2004, and is under construction at 
the time of this writing. 



4. COBRA (Cadmium telluride O-neutrino Beta decay Research Apparatus): 



COBRA (jZuberl . I2OOII ) is an R&D program developing CdZnTe semiconductor crystal 
detectors (usually referred to as CZT detectors) for f3(3{0v). CZT detectors contain 9 isotopes 
that undergo (3(3. The isotopes "*^Cd, ^^OTe, ^^^Cd, ^°Zn, and ^^STe all undergo (3- (3- , 
whereas the isotopes ^^Zn, ^o^Cd, ^^^Cd, and ^^oxe undergo /3+EC, or EC-EC. The 
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critical nucleus is ^^^Cd as it has the highest Qf^p and therefore the greatest sensitivity to 
/9/9(0i^). Although the Qj3p of ^^°Te is also high, it is below that of ^^^Cd and therefore the 
Cd (3(3{2i') would create a significant background for any /5/?(0i/) study of Te. The isotopic 
abundance of ^^^Cd is 7.5%, however, necessitating enrichment for a competitive experiment. 
The existence of a number of /?+EC isotopes within the detector and the high granularity 
of the experiment permit a variety of studies in addition to the primary /?/5(0i/). 

The proposal is to operate 64,000 1-cm^ CZT detectors (^^^Cdo.gZno.iTe) with a total 
mass of 418 kg. The detectors would be fabricated with Cd enriched to 90% in isotope 116, 
so about 44%, or 183 kg, of the detector mass is the isotope of interest. The energy resolution 
(FWHM) of the co-planar grid (CPG) detectors is better than 2% at the endpoint. If the 
background can be kept below 1 count/ (keV t y), a half life sensitivity of better than 10^^ 
y is anticipated. To improve the resolution of the detectors, cooling and new grid designs 
are being investigated. The collaboration has operated a 2 x 2 detector prototype and is 
assembling a 4 x 4 x 4 prototype array at Gran Sasso (Fig. [T5|) . with the first 16 detectors 
running since s ummer 2006 (Fig. [T6ll . 

Initial data (jGo/31ing et al\ . 120051 ) for the 2x2 detector array, including a measurement 
of the 4-fold forbidden 13 decay of ^^^Cd and new limits on positron decay modes of (3(3, have 
been published. The collaboration is also studying pixellisation of the detectors to improve 
spatial resolution within a single detector for background rejection by particle identification. 



5. MOON (The Molybdenum Observatory of Neutrinos) 



The MOON experiment is unique because it has two applications. First, it is a tracking 
detector to search for the double-beta decay of ^°'^Mo and other isotopes. The energy and 
angular correlations of individual (3 rays are used to identify the z/-mass term of the /3/3(0z/) 
decay. In addition it will detect solar neutrinos, from the low-energy pp-neutrinos all the 
way up to ®B neutrinos. This du al fun c tion results from t he un ique nuclear structure of the 



triplet i""Mo, ^""Tc, and ^""Ru (lEiiril 



2005 



Ejiri et al. 



concen trate on MOON as a double-beta decay experiment (lEjiril . 



20001) . In this discussion we will 



2006b 



Nakamura et al. 



2006a 



a) 



The MOON collaboration involves members from 12 institutions in Japan, the Czech 
Republic, Russia and the US. The goal is to achieve an effective sensitivity to (m^^) of ~50 
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meV with one ton of ^°°Mo. This isotope has a high Qpp (3034 keV), above the highest 
energy and troublesome 2615-keV 7 ray in ^°^T1. It has a correspondingly large phase-space 



factor, about 7 times larger than that of Ge, and a natural isotopic abundance of 9.( 



70. 



While its QRPA nuclear matrix element is slightly smaller than that of ''^Ge, its predicted 



rate is 3.6 times faster ( Rodin et al. 



20061 ). On the other hand, for the value of Qpp used 



in the calculations of Civ itarese and Suhonen , this ratio is 5.84. Neither the old nor recent 



Shell Model calculations (ICaurier et al. 



1996 



Poves et al. 



20061 ) treat the case of ^°°Mo 



The above two recent QRPA calculations differ by a factor of 1.6, and with other values of 
Qpp, as much as a factor of 4. This makes a great difference in the predicted effectiveness of 
^°°Mo [3(3 experiments. Thus depending on the chosen value for Qpp, the predicted l3[3{^u) 
rate might vary by a factor of 4. It should also be pointed out that the MOON technique, 
where the source does not comprise the detector, could also be used with other isotopes, 
^^Se and ^^°Nd, for example. The choice is made by considering Mqi,, Q/3/3, the /3/3(2z/) rate, 
and other characteristics of a given isotope. 



The experimental technique is based on the ELEGANT V design (lEjiri et ali ll991al Jbl). 
The apparatus consists of multi-layer modules, each with a thin film of enriched Mo, sand- 
wiched between two posit ion- sensitive detector planes (thin MWPC chambers or thin PL 
fibers) and two solid scintillator (plastic) plates, with all other modules acting as an active 
shield. Precise localization of the two /3-ray tracks enables one to select true signals and reject 
background. In fact, background from radioactive impurities and from neutrons associated 
with cosmic-ray muons is evaluated by Monte-Carlo simulations and found to be acceptable 
(< 0.3 /(t y)). One of the key challenges in this type of experiment is achieving adequate 
energy resolution to be able to recognize (3(3{Qv) events from background, and also from the 
irreducible background from (3(3{2v) events. The half life for (3(3{2v) of ^°°Mo was measured 
ve ry accurately by the NEMO experiment as Tlj^ = 7.11 ± 0.02{stat.) ± 0.54(s?/st.) x 10^^ 
y (lArnold et all l2005al ). The theoreti c ally p redicted of ^°°Mo, for a {mpj^) = 50 meV, 
is as long as 3.8 x 10^^ y ( iRodin et a/.l . l2006l ). although other Mq^ calculations indicate 
nearer to 10^^ y. The (3l3{2v) rate is then between 7 and 8 orders of magnitude faster than the 
predicted (3(3{Qv) rate. Obviously, energy resolution is a key issue in all s uch experiments 



In the case of MOON, it is being addressed by an intensive R&D program (INakamura et al. 



2006ah . 



The MOON-I prototype has the general features shown in Fig. [T71 The three 48 cm x 
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48 cm molybdenum films are about 47 g each and are enriched to 94.5% in ^"^Mo. They are 
supported and covered by aluminized Mylar films, 6 /xg/cm^. This prevents optical cross talk 
between the scintillators. The six plastic scintillators are 53 cm x 53 cm x 1 cm, and are 
viewed by 56 Hamamatsu, R6236-01 KMOD photomultiplier tubes that are relatively low 
in ^°K. They are optically coupled to the plastic scintillators with silicon cookies, which are 
made from silicon rubber 3 mm thick. The general geometry of the MOON-I configuration 
is shown in Fig. [181 The pattern of PMT coupling to the scintillators shown in Fig. [T71 
allows the hit- scintillator to be identified by the PMT hit pattern. The measured energy 
resolution is a = 4.8 ± 0.2%(AE/E) at 1 MeV. This would yield a (x ^ 2.9% and FWHM of 
approximately 200 keV at the Qp^^ of 3034 keV. The collaboration is attempting to improve 
it to cr ~ 2.2% (as they have done in a small prototype) by adjusting the position dependent 
response. 

The MOON-I experiment was performed with the detector inside the EL EGANTS-V 



shield in the Oto underground Labotatory with an overburden of 1300 m.w.e (lEjiri et al. 



1991a 



20011 1. At this level the muon flux was measured to be 4 x 10 ''yu/ (cm^ s). There are 



Nal(Tl) detectors above and below the MOON-I array to detect background 7 rays. The 
shield consists of 15 cm of lead outside of 10 cm of oxygen-free, high conductivity copper. 
The shield was flushed with dry N2 that reduced the activity from Radon to a level of 125 
mBq/m^. 

The MOON research and development program continues with the main goals of reducing 
background as indicated by simulations (which have not yet been reported), improving the 
energy resolution, and improving the position resolution for tracking. The sensitivity of any 
^o^Mo experiment to {m/sp) cannot really be determined until more reliable nuclear matrix 
elements are available. Experimental studies with charge exchange r eactions are in progress 
for s^Se, i°0Mo, and ^^"Nd, which are candidates for use in MOON jEiiri booof l. 



6. NEMO-3 and SuperNEMO 



NEMO-3: The SuperNEMO detector flBarabashl . l2004f ) wil l be a trackiiig detect or, as 



is its currently operating predece ssor, the NEMO-3 detector f)Arnold et al\ . l2005bl ). We 



discuss the NEMO-3 experiment (lArnold et al 



2005a 



Barabash 



2006bD and its results to 



provide background. NEMO-3 is the third generation of NEMO j3(3 detectors. An artist's 
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conception is shown in Fig. [191 The detector has 20 segments of thin source planes, with a 
total area of 20 m^, that can support about 10 kg of source material. It has a 3-dimensional 
readout wire drift chamber, with 6180 cells that operate in the Geiger mode for tracking. The 
detector gas is He, with 4% ethyl alcohol, 1% argon, and 0.1% H2O. The tracking volumes 
are surrounded by 1940 plastic scintillator-block calorimeters. The scntillator detectors 
operate with thresholds of ~30 keV, and have efficiencies of 50% for 1-MeV 7 rays. The 
energy resolutions range from 11% to 14.5% FWHM at about 1 MeV. The resolution is one 
parameter that will be given a lot of attention in the design and construction of SuperNEMO. 
Good energy resolution is critical for the discovery of processes with long half lives. The 
source planes hang vertically in a cylindrical geometry inside of a magnetic solenoid that 
generates a ~25 G magnetic induction field. Tracking in the magnetic field allows the 
differentiation between electron and positron tracks, with only a 3% chance of confusing the 
two. 

The detector is surrounded by 18 cm of low-background iron to reduce the external 7- 
ray fiux. Fast neutrons from the laboratory environment are suppressed by an external 
shield of water, and by wood and polyethylene plates. NEMO-3 is located in the Modane 
Underground Laboratory in Frejus, France with an overburden of 4800 m.w.e.. The air in 
the experimental area is constantly fiushed and has a radon-free purification system serving 
the detector volume. It reduces the radon to ^200 mBq/m^, and has a capacity of 150 
m^/h. The measured radon activity in the detector was 4-5 mBq/m^. 

One of the great advantages of the NEMO-3 detector is that it can make measurements 
on many different nuclei, and in fact has made measurements on seven isotopes, and two 
blanks, simultaneously. The isotopes were ^°°Mo (6.914 kg), ^^Se (932 g), "'^'Cu (621 g), 
"'^'Te (491 g), i30Te (454 g), ^^^cd (405 g), ^^^Nd (37 g), ^^Zi (9.4 g), and ^^Ca (7 g). The 
"°*Te and ""*Cu were used to measure the external background. 

The measur ement of the decay of ^°°Mo made with NEMO-3 truly sets a standard for such 



measurements flArnold et fli.l . l2005al ). Th e result wa s : TfY^ = 7 .ll±0.Q2{stat)±QM{syst) x 
10^^ y. New bounds were set on /?/3(0z/) (IBarabashl . l2006bl ): T^^^ > 5.8 x 10^^ y for i°°Mo, 
and 2.1 x 10^'^ y for ^^Se (90% confidence level). The bounds on the effective Majorana 
mass of the electron neutrino are {mpj^) <0.6-2.7 eV for ^'^''Mo and 1.2-3.2 eV for '^^Se. The 
NEMO-3 detector has been operating since February 2003. The collaboration projects that 
the sensitivity to T^J^ of ^°°Mo will be 2 x 10^"^ y by the end of 2009. This would reduce 
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the bound on {mpp) to <0. 34-1. 34 eV. It is clear that by that time, the detector will have 
been operating for seven years and the resulting sensitivity will not be competitive with 
next-generation experiments that are designed to probe the inverted-hierarchy mass scale. 
To remedy that situation, the SuperNEMO Collaboration was formed. 

SuperNEMO: This proposed experiment is a vastly expanded tracking chamber of a 
modular design. The parent isotope will either be ^^Se or ^^°Nd. At this time there does 
not exist a source of kilogram quantities of Nd isotopically enriched in ^^°Nd, whereas ^^Se 
has been enriched by the gas centrifuge technique in Russia. In Table |V] below we briefly 
compare the experimental parameters of NEMO-3 to those of SuperNEMO. 

At the time of this writing, SuperNEMO is a collaboration of 26 institutes in 11 countries. 
An artist's conception of the detector modules is shown in Fig. [201 The detector will comprise 
20 modules. Each module is designed to hold 5 kg of enriched isotope; each has 12 m^ of 
tracking volume, with 3000 channels of readout and 1000 photomultiplier tubes. The total 
detector will have 60,000 channels of drift chamber readout, and 20,000 photomultiplier 
tubes. 

One of the key issues is the background from the decay of ^^^Bi to ^^^Po in the ^^^^U 
chain and ^^^Bi to ^^^Po in the ^ss^j]^ chain. These are referred to as the BiPo backgrounds 
and a prototype-test detector named BiPo is installed in the Laboratorio Subteranio de 
Canfranc (LSC) in Canfranc Spain to test the source foils. The BiPo detector comprises two 
modules. The intent is to test each source foil for the BiPo backgrounds before installation 
in SuperNEMO. The plan is to measure a 5 kg foil every month with a sensitivity of ^os^^ 
< 2/iBq/kg and ^^^Bi < lO/xBq/kg. 

The proposed schedule is as follows: BiPol will run until 2008; BiPo will be under 
construction from mid 2007 until mid 2009. The final BiPo will be installed in LSC at the 
beginning of 2009. The construction of the first SuperNEMO module will start in early 2009 
and be completed in mid 2010. The BiPo tests will continue throughout the construction of 
SuperNEMO. The full detector is scheduled to start operating at the beginning of 2012. 



59 



7. CANDLES ( CAlcium fluoride for studies of Neutrinos and Dark matters by Low Energy 
Spectrometer) 



The proposed CANDLES experiment (lUmehara et all l2006l ) is based on un-doped CaF2 
scintillators surrounded by a 4tt liquid scintillator live shield. The scintillation light from 
pure CaF2 crystals is in the ultra-violet region of the spectrum and has a quite long decay 
time (1 /is) although the liquid scintillator has a very short decay time (~10 ns). One can 
employ this difference to reject background signals from the liquid scintillator. A wavelength 
shifter will be added to the liquid scintillator to match the response of the photomultiplier 
tubes. The CaF2 crystals are made with natural- abundance calcium which is 0.187% ^^Ca. 
The endpoint energy of ^^Ca, 4.27 MeV, is the largest of all the PP candidates, well above 
the highest 7 ray from natural radioactivity, namely the 2615-keV 7 ray in the decay of 
208'jj_ jg interesting to note that ^^Ca is a good candidate for shell-model analysis of the 
f](3 matrix elements. 

The CANDLES project follows the ELEGANT VI experiments performed with CaF2 
scintillation detectors made with europium-doped CaF2(Eu). In the ELEGANT VI detector 
there were 23 CaF2(Eu) crystals, each 4.5 x 4.5 x 4.5 cm'^ (290 g), surrounded by 46 CaF2 
(pure) crystals that acted as a shield and as light guides. There was 6.7 kg of CaF2 containing 
6.4 g of ^^Ca. This in turn was surrounded by 38 CsI(Tl) crystals, each 6.5 x 6.5 x 25 cm^, 
to act as a live shield. The entire array was placed inside a bulk shield of copper, an air-tight 
box, lead, Li-loaded paraffin and Cd, and B-loaded water. The experiment was located in 
the Oto Cosmo Observatory with an overburden of 1400 m.w.e.. The energy resolution was 
3.1% at 4271 keV, and the detection e fficiency was fti49%. The experiment resulted in a new 



bound for T ^U^^Ca) > 



calculations (jPoves et al 



. 8 X 1 0^^ y (lOgawa et a/.l . l2004l ). According to recent shell model 



20061 ). {mpji) < 23 eV. It will require several ton years for such 



an experiment to be competitive. Currently CANDLES-III is under construction with 60 
crystals for a total mass of 191 kg at a sea level lab in Osaka (see Fig. [2T1) . 

The proposed CANDLES detector will comprise several tons of CaF2(pure) detectors, 
with almost no background, immersed in liquid scintillator that will act as a veto detector. 
The target parameters for a 6.4-ton CANDLES are < 3 /iBq/kg of both 212,2143^^ 3 5% 
FWHM resolution, and 6 years of counting with a total background of 0.3 per year under 
the expected /3/3(Oz/) peak. The target sensitivity for this phase is {m^p) ^0.1 eV. According 
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to shell model calculations (jPoves et al. 

X 10^6 y. 



20061 ). this would require a sensitivity of T^j2 ~ '''•2 



The shell model calculations predict a value for the nuclear structure factor F/v (^^Ca) 
= 3.61 X 10-^7 y. The same model predicts Fjv(™Ge) = 3.22 x lO^^^/y and Fjv(^^°Te) = 
2.13 X 10^^^/y. One can immediately appreciate the challenge of performing an experiment 
with natural- abundance calcium. One would need about 400 times as many Ca atoms as 
Ge atoms, enriched to 86% in ^^Ge, to have the same decay rate. In comparison with ^^°Te, 
one would need about 1000 times the natural abundance Ca atoms as natural abundance 
of Te atoms to have the same theoretical decay rates. The collaboration is investigating 
techniques for enriching Ca. 



8. Other Proposals 



CARVEL (jZdesenko et alV l2005l ) (CAlcium Research for VEry Low neutrino mass) is 
a proposal to use isotopically enriched ■*^CaW04 crystal scintillators to search for /5/5(Oz/) 
of ^^Ca. These scintillators, with natural calcium, have been tested for energy resolution, 
internal background and pulse-shape discrimination that enables the discrimination between 
alpha and gamma ray events. The energy resolution at 2615 keV is 3.8% (~100 keV) which 
compares well with tracking detectors. The total background between 3.8 and 5.4 MeV, 
after the elimination of a background by pulse-shape discrimination, was 0.07 counts/(keV 
kg y). With these parameters it was determined that with ^100 kg of ^^CaW04, enriched 



to 80% in ^Ca, the half-life sensitivity would be T^^g 



Ca) ~ 10^^ y. T 



l is cor responds to 



2006|). With one 



a sensitivity in effective neutrino mass of {mpp) ~ 55 meV (iRodin et ali 
ton of this material, it was estimated that a half-life sensitivity of 10^^ y could be reached, 
corresponding to {m^p) ~ 0.02 eV. The difficulty with this proposal is that there is no source 
of kg quantities of Ca enriched in '^^Ca. While in principle Ca can be enriched by Atomic 
Vapor Laser Isotope Separation (AVLIS), the effort and cost scale as the ratio of the final 
to initial abundance. This cost would be more than 35 times that for enriching an element 
with a 7% natural abundance [ e.g. ^^Ge). This is a challenge for a ton-scale experiment. 

DCBA ( jlshihara et ali . i200Q) (Drift Chamber Beta-ray Analyzer) is designed to search 
for the /5/3(Oz/) of ^^°Nd. It comprises tracking chambers, a solenoid magnet and a cosmic-ray 
veto detector. The source foils are vertical with a series of high- voltage anode wires parallel 
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to the source and the magnetic field. As the /3-particles are emitted they are given an extra 
velocity perpendicular to the magnetic field, and are curved as they pass through the 1-atm 
He gas mixture. The anode wires detect the drift electrons, while the cathode wires across 
the chamber gather the drifted ions. The detector can measure the momentum of each (3 
particle and can determine the vertex. Drift electrons generate avalanche electrons on the 
anode wire plane. The pattern of avalanche ions and the time of each signal are recorded 
with fiash ADCs. From these data, the x-position of the origin of the track is determined. 
The y-position is fixed by the responding anode wire. The position on the anode wire of 
the avalanche determines the z-position with pickup-wires, which are located transversely 
near anode wires and pick up signals induced by the avalanche. Each electron track is 
reconstructed in 3-dimensions, and the momentum and kinetic energy are obtained from the 
curvature and pitch angle. DCBA-I will contain natural Nd and will have ~ 3.8 x 10^^ atoms 
of ^^°Nd (5.6% natural abundance). It is expected to have a background of < 0.1 counts/(keV 
y). It will have an energy resolution of :^100 keV. DCBA-II is planned to have 4.5 x 10^^ 
atoms of ^^°Nd, in the form of Nd203 with Nd enriched to 90% in ^^° Nd. The expected 
half-life sensitivity is 4 x 10^^ y. According to recent QRPA calculations ( IRodin et all l2006l ) 
this would correspond to {mj3p)r^ 0.12 eV. The R&D for this detector is being carried out at 
KEK with a test prototype DCBA-T. On the basis of the that R&D a new pro j ect h as been 
proposed, temporally named MTD (Magnetic Tracking Detector) (jlshiharal . 120071 ). The 
detection principle is the same as in DCBA, but the amount of source installed is much 
larger. One module of MTD will contain natural Nd and will have ~ 6.7 x 10^^ atoms of 
^^°Nd (5.6% natural abundance) in the first experimental phase. The natural source plates 
will be replaced with enriched material in the second phase. The enrichment of ^^°Nd is 
planned on a large scale w ith an internationa l collaboration. 

I d inst all 1000 kg of "^CdW04 crystals in 



The CAMEO proposal (iBeUini et al 



2001) 



wou. 



the liquid scintillator of BOREXINO flBeUini et al 



19961 ) after its solar neutrino program is 



complete. The described program would begin (CAMEO-I) with thin (^15 mg/cm^) metal 
isotopically-enriched Mo sheets sandwiched between plastic scintillator and installed within 
the scintillator of the BOREXINO Counting Test Facility (CTF). Better sensitivity could be 
achieved, however, with an active source configuration of CdW04 crystals. For CAMEO-II, 
100 kg of crystals would be installed in the CTF. Measurements of the intrinsic background 



of the crystals lead to an estimated sensitivity of ~ 10 y. Installing 1 ton of crystals 
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into the BOREXINO detector could allow a sensitivity of 10^^ y. Because of the focus on 
solar neutrinos at BOREXINO, this program is not currently active. 

Another proposal th at would exploit the BOREX I NO d etector or the CTF is to dissolve 



20011 ). In the most ambitious design, a 



^^^Xe in the scintillator (jCaccianiga and Giammarchil . 
vessel of radius 2.7-m containing Xe and scintillator is itself contained with the BOREXINO 
vessel (radius 4.25 m) filled only with scintillator. Outside this scintillating region is a 
pseudocumene buffer. About 2%-by-weight-Xe is soluble in scintillator with an appreciable 
temperature dependence, so an estimated 1565 kg of Xe could be used at 15° C. Simulations 
indicate that phototube activity and /3/3(2z/) will dominate the background, resulting in a 
T^^2 sensitivity of 2 x 10^^ y. This program i s currently inactive. 

The XMASS experiment ( iTakeuchil . |200J) (Xenon detector for weakly interacting MAS- 
Sive particles) is a dark-matter search using liquid Xe as a target for WIMP detection via 
nuclear scattering. The present detector uses 100 kg of natural Xe with a fiducial volume 
of 3 kg viewed by phototubes. It is planned to expand this to 1 ton and even possibly to 
20 tons. This final configuration would have a fiducial volume of 10 tons and would be used 
to study dark matter and solar neutrinos. With such a large mass of Xe, the sensitivity 
to /5/5(Oz/) is clearly of interest. The self shielding of the Xe is highly effective at reducing 
the backgrounds at low energies (<500 keV) of interest to solar neutrinos and dark matter. 
However, it is less successful at the higher energy appropriate for (5(5{^v). Present efforts 
are aimed at reducing the backgrounds for the dark matter search. Results from this 100- 
kg prototype have shown that the dark matter configuration will not work for which 
will require a dedicated design. Work toward such a design is not a current priority of the 
collaboration. 

There are two proposals to use Ce-doped GSO scintillating crystals (Gd2Si05:Ce) to 
search for /3/3(Oi/) in ^^°Gd [Qnn = 1-73 MeV). In one research effort, a 1.744-kg GSO 
crystal was used (jWang et al\ . l2002l ) to study backgrounds. This crystal, coupled to a 
2-inch phototube, was wrapped in Teflon and black vinyl tape. The whole crystal was 
inserted into a well-type Nal crystal. The system was contained within a Pb shield and a 
plastic-scintillator cosmic-ray veto. The whole setup was situated in the bottom floor of 
a 7-story building (12 m.w.e. overburden). The crystal was found to have a substantial 
contamination from U- and Th-chain isotopes, which were not in equilibrium. In addition, 
a decays of ^^^Gd were present. This program is no longer being pursued. GSO crystals 
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were also used in an experiment at the Solotvina Underground Laboratory (IDanevichl . l200ll ) 
at a depth of ^1000 m.w.e.. In this test, a 635-g crystal was joined to a phototube by a 
18.2-cm plastic light guide. A passive shield of Cu, Hg and Pb surrounded the detector. 
Internal radioactive contaminants within the crystal were the limitation of this experiment, 
although the detector was cleaner than that discussed above. This collaboration proposed 
putting GSO crystals within a large volume of liquid scintillator (an idea similar to that of 
the CAMEO project) to help control the background. Background levels will be a challenge 
for this technology. Although GSO crystals are costly, the high isotopic abundance of ^^°Gd 
(22%) would permit the use of natural Gd, which would be a great relative cost reduction 
compared to other proposals, assuming that pulse shape techniques could be used to suppress 
backgrounds. A modest R&D progr am continue s on this project. 



Final 



y, there is a l so a p roposal (IChenl . |2005| ) to use the Sudbury Neutrino Observatory 



(SNO) (IBoger et all 120001 ) for PP, now that it has completed its solar neutrino studies. 



The plan for the SNO detector is to replace the heavy water with scintillator, beginning a 
program of low-energy solar-neutrino studies and long-baseline-oscillation studies of reactor 
neutrinos, geoneutrinos, and PP. The parent isotope, situated in the liquid scintillator 
either as nanoparticles, dissolved chemicals, or as absorbed gas, can be used as a source, 
with a Nd-loaded liquid scintillator as the leading candidate. The scintillator configuration 
is referred to as SNO+. A 1% loading of natural Nd in the scintillator would provide 
10 tons of Nd, corresponding to more than 500 kg of ^^°Nd isotope. Simulations of this 
configuration indicate that PP{2v) will be the dominant background because of the limited 
energy resolution. However, if {mpp) is in the degenerate mass-scale region, the extremely 
high count rate from this large isotopic mass would permit a statistical separation of PP{2v) 
and PP{Qv). A development program for SNO+ is proceeding. 

There are a number of research and development projects involving other crys- 
tals that could be candidate s for PP detectors. The cryst als under i nvestigation are 



CdWp4 fiBardeUi et al 



200fih. PbW O. f 



also being considered (INygrenl . 



Danevich et al. 



2005aD, and ZnW04 (iDanevich et all [2005^ A detector based on high pressure Xe gas is 



20o"fil ). YAG:Nd fiDanevich et all 



20071 ). 
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VIII. CONCLUSIONS 



The future for (3(3 is exciting. Technologies in hand will allow us to measure at intrigu- 
ing neutrino-mass levels. Technological progress is rapid, and developments in theory are 
improving our ability to interpret the measurements. Strong (3(3 programs in both particle 
and nuclear physics are the result. 

The answer to the question: "/s the neutrino its own anti-particle?^ , is critical for theories 
of particle mass. It is also needed to help uncover the reasons matter dominates over anti- 
matter in our Universe. The neutrino mass will not only tell us about the high scale at which 
the Standard Model breaks down {e.g., through the see-saw), but also will have imphcations 
for the large-scale structure of the Universe. Finally, lepton-number violaton is significant 
in its own right. 

(3(3{0i') is the only practical method for investigating the particle/anti-particle question. 
And if the neutrino is its own anti-particle, (3(3{0i') will have the best sensitivity to neutrino 
mass of any laboratory technique. The (3P program outlined here, consisting of several 
experimental results, will therefore greatly influence the wider-particle physics endeavor. 

Questions about how best to calculate the nuclear matrix elements have led to develop- 
ments in nuclear theory, despite the preoccupation of most theorists with other problems. 
Large-scale shell-model codes are increasing in power, in part because of the need for better 
matrix elements. And a variety of experiments on nuclei are planned to constrain those 
same matrix elements. We should see some reduction in their uncertainties. 

Experimental progress in (3(3{Qv) has also led to improved material-purification techniques 
such as Cu electroforming, and to ongoing improvement in assay techniques such as direct 
7-ray counting and inductively-coupled plasma-mass spectroscopy. The need for effective 
detectors has led to many improvements in semiconductor {e.g. Ge, CZT) and bolometer 
technologies that are now finding application in basic science, medicine, and homeland se- 
curity. As researchers strive for 1-ton (3(3 experiments, advances in areas such as isotope 
enrichment and detector production should follow. 

Overall, the program in (3(3 has been and will continue to be very fruitful. 
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FIG. 1 A generic level diagram for (3f3. 





Vm X 




FIG. 2 Feynman Diagrams for (3(3{2v) (top) and (3(3{^v) (bottom). 




FIG. 3 Majorana propagator resulting from (3(3{0i^) amplitude (jSchechter and Valle 



1982 ). 
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FIG. 4 The effective Majorana mass {mpp) as a function of the mass of the hghtest neutrino, 
miightest- In making the plot, we have used the best fit values for the parameters in Table HI The 
filled areas represent the range possible because of the Majorana phases and are irreducible. If 



one in corporates the uncertainties in the mixing parameters, the regions widen. See 



Bilenky et al 



( 20041 ) for an example of how the mixing parameter uncertainty affects the regions. 




FIG. 5 The distribution of the sum of electron energies for (3(3{2v) (dotted curve) and (3(3{Qv) (solid 
curve). The curves were drawn assuming that Tqi, is 1% of and for a l-cr energy resolution of 



70. 
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FIG. 6 The distribution of the electron energies taken one at a time for f3(3{Qv). The soHd curve is 
for light-neutrino excha nge and the dott ed curve is for right-handed currents in both the leptonic 



and hadronic currents. (Doi et al 



1983 ). 




FIG. 7 The CUORICINO structure. Figure courtesy of the CUORE collaboration. 
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FIG. 8 The CUORE structure. Figure courtesy of the CUORE collaboration. 
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FIG. 9 The atomic level structure of a Ba ion indicating the levels of interest and a picture of the 
emission of a lone ion within an ion trap. Figure courtesy of the EXO collaboration. 
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FIG. 10 An artists conception of the EXO-200 experiment. Figure courtesy of the EXO collabo- 
ration. 




FIG. 11 The left panel shows a concept drawing of 1 of 19 strings of Ge detectors that would be 
contained within a Majorana module shown in the right panel. Figure courtesy of the Majorana 
Collaboration. 
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FIG. 12 The top panel shows the waveform of a candidate single-site energy deposit in a Ge 
detector. The middle panel shows a candidate multiple-site energy deposit of similar energy. The 
bottom panel shows how the spectrum in the region near 1.6 MeV changes as pulse shape analysis 
and segmentation cuts are employed. Figure courtesy of the Majorana Collaboration. 
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GERDA design 




FIG. 13 An artist's conception of the GERDA experiment. 
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FIG. 14 Artist's conception of the GERDA Phase II tower. 




FIG. 15 An artist's conception of COBRA 64-detector setup. Figure courtesy of the COBRA 
collaboration. 




FIG. 16 A photograph of a detector holder containing 16 detectors, courtesy of the COBRA 
collaboration. 




FIG. 17 A schematic of the MOON-I scintillator layout, courtesy of the MOON collaboration. 
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FIG. 18 A schematic of the MOON-I configuration, courtesy of the MOON collaboration. 
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FIG. 20 A schematic of a SuperNEMO module. Figure courtesy of the NEMO collaboration. 




FIG. 21 A photograph of the CANDLES-III installation. Figure courtesy of the CANDLES col- 
laboration. 



Tables 



FAB 



lE I Neutrino mixing parameters as summarized by the Particle Data Book ()Yao et al. 



(1200a)) based on the individual experimental reference reporting. The limit on {mp) and E are 



based on the references given. The {mpp) limit comes from the Ge experiments. Th e parameter 



value s would be slightly different if determined by a global fit to all oscillation data (jFogli et al. 



20061 ) 



Parameter 



Value 



Confidence Level 



Reference 



sin2(26l23) > 0.92 

sm2(2(9i3) < 0.19 

Amii S.Oi^l X 10-5 



-0.3 







lAmial 2.4t°| x 10-^ eV^ 
{mp) <2eV 
{mpp) < 0.7 eV 



70 

95% 



Aharmin et al. (2005 ) 

Ashie et al. (2005) 
Apollonio et al. (1999) 
Aharmin et al. (2005) 
Ashie et al. (2004) 
Kraus et al. (2005 ): Lobashev et al. (1999) 
Aalseth et al. f 2002a ): Klapdor-Kleingrothaus et al. (2001a) 
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S <2eV I m% — 

"Using the matrix element of lRodin et al\ (j2006f ) 



Elgarov and Lahov (2003 ) 



TABLE I I A lis t of the values of ^'^'^ various isotopes. These values are recommended by 
BarabashI (2006a) a s the best interpretation of the experimental data. One should heed discussion 



m 



Barabash ( 



2006al ) before using the values in a quantitative way. 



Isotope 


Tf- (years) 


Isotope 


Tf- (years) 


48Ca 


(4.2t2-i) X 1019 


i28Te 


(2.5 ±0.3) X 10^4 


76Ge 


(1.5 ±0.1) X 10^1 


i30Ba EC-EC(2zy) 


(2.2 ±0.5) X 10^1 


82Se 


(0.92 ± 0.07) X 10^0 


isoTe 


(0.9 ±0.1) X 10^1 


96 Zr 


(2.0 ±0.3) X 10^9 


isoNd 


(7.8 ±0.7) X 10^8 


iooMo 


(7.1 ±0.4) X 10^8 


238 u 


(2.0 ±0.6) X 10^1 


ii^Cd 


(3.0 ± 0.2) X 10^9 







TABLE III A list of recent f3f3{Qv) experiments and their 90%-confidence-level (except as noted) 
limits on Ty^- The {mpp) limits, if provided, are those quoted by the authors, who each made 
choices about which calculated Mq^ to use. 



Isotope 


Technique 


^1/2 


(m^^) (eV) 


Reference 


48 Ca 


CaF2 scint. crystals 


> 1.4 X 10^2 y 


< 7. 2-44. 7 


Oeawa et al. (2004) 


76Ge 


^"'^Ge det. 


> 1.9 X 10^5 y 


< 0.35 


KlaDdor-Kleinerothaus et al 


76Ge 


^"'■Ge det. 


(2.23l°:44) X 1025 y [la) 

> 1.57 X 1025 y 


0.32 ±0.03 


KlaDdor-Kleinerothaus and Krivc 


76Ge 


^"'■Ge det. 


< (0.33 - 1.35) 


Aalseth et al. f 2002a 


82Se 


Thin metal foils and tracking 


> 2.1 X 1023 y 


<(1.2-3.2) 


Barabash ('2006bl 


iooMo 


Thin metal foils and tracking 


> 5.8 X 1023 y 


<(0.6-2.7) 


Barabash ('2006bl 


"6Cd 


ii6CdW04 scint. crystals 


> 1.7 X 1023 y 


<1.7 


Danevich et al. (200c 


i28Te 


geochemical 


> 7.7 X 1024 y 


<(1.1-1.5) 


Bernatowicz et al. (19 


i30Te 


Te02 bolometers 


> 3.0 X 1024 y 


<(0.41-0.98) 


Arnaboldi et al. (200 


i36Xe 


Liq. Xe scint. 


> 4.5 X 1023 y-^ 


<(0.8-5.6) 


Bernabei et al. (2002 




Thin metal foils and tracking 


> 3.6 X 1021 y 




Barabash ^2005^ 
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"See footnote 4 in iBernabei et al\ ()2002D 
TABLE IV A summary list of the /?/3(0z/) proposals and experiments. 



Experiment 


Isotope 


Mass 


Technique 


Present Status 


Reference 


P A 1\TDT TPQ 


>^a 


few tons 


CaF2 scint. crystals 


Prototype 


Umehara et al. (1 




Ua 


1 ton 


CaW04 scint. crystals 


Development 


Zdesenko et al. (2 




116pj 


^lo Kg, 


v^zj ± semicuiiQ. clcl. 


r lototype 


Zjuuei \ zuui I 


PT TPT3 TPTMP 


le 


4U. ( Kg 


Te02 bolometers 


Running 


Arnaboldi et al. (1 


CUORE 


isoxe 


741 kg 


Te02 bolometers 


Proposal 


Ardito et al. (20 


RlPT5 A 


150i\T„ 

i\e 


2U kg 


^'^^Nd foils and tracking 


Development 


Ishihara et al. (2{ 




136Xp 
-A.e 


ZUU Kg 


ijiC[. vVe ± Jr \^ / SClllL. 


V.^Oilbtl UCXlOIi 


r^ieuKe I ZUU / j 




136Yp 


1 1 n f 


Ijiq. JVG Ir L//scmt. 


Proposal 


uaniiov ei at. ( zuu 


PFM 


76 Pp 


1 ton 


vjc U.t:l. Ill llt^. lllLlUgcll 


lilciCLl VtJ 


ZjU-tioCllKO C.fc tit. I Z 




vje 


1 ton 


KjKi Qei. Ill ncj. nitiogen 


iiiacLive 


-LviauQoi-jA.iemiiioinaus e 






--SoD Kg 


^'^'"Ge semicond. det. 


Construction 


Schonert et al. (2 


PQP 


160 p J 


z ton 


Gd2Si05:Ce crys. scint. in liq. scint. 


Development 


JJanevicn et al. ( zuuu ) . wan 


Majorana 


76Ge 


120 kg 


^"■^Ge semicond. det. 


Proposal 


Gaitskell et al. ( 2 


MOON 


iooMo 


1 t 


'^"'"Mo foils/scint. 


Proposal 


Nakamura et al. (2 


SNO++ 


150Nd 


10 t 


Nd loaded liq. scint. 


Proposal 


Chen (2005) 


SuperNEMO 


82Se 


100 kg 


enrgg fgils/tracking 


Proposal 


Barabash (200^ 


Xe 


i36Xe 


1.56 t 


^"'^Xe in liq. scint. 


Development 


Caccianiea and Giamma 


XMASS 




10 ton 


liquid Xe 


Prototype 


Takeuchi (2004 


HPXe 


i36Xe 


tons 


High Pressure Xe gas 


Development 


Nveren f2007' 



TABLE V A comparison of the SuperNEMO design parameters with those of the NEMO detector. 



NEMO 


Detector /Experiment 


SuperNEMO 


ioomo 


Isotope 


82sg 150^^ 
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7 kg 



o/o 

< 20 ^Bq/kg 
^i^Bi < 300 /xBq/kg 
8% at 3 MeV 
2 X 10^4 y 
0.3-1.3 eV 



Source Mass 
(3(3{0iy) Detection Efficiency 
External Background 
In the source foil 
Energy Resolution, FWHM 
T^J^ Sensitivity 
{f^jSp) Sensitivity 



100-200 kg 

^5^30% 
208X1 < 2 /xBq/kg 
For »2Se:2i4Bi < 10 /xBq/kg 
4% at 3 MeV 

> 10^6 y 
<0.05-0.1 eV 
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